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SCIENTIFIC  OBJECTIVE: 

The  objective  of  this  work  is  to  use  mostly  optical  techniques  to  study  quantum  tunneling  struc¬ 
tures,  quantum  wells,  semiconductor  super  lattices  and  alloy*  grown  by  MBE  (and  occasionally  other 
techniques  such  as  multitarget  i£.  sputtering).  Of  primary  interest  are  the  effect  of  super  lattice  layer¬ 
ing  and  alloy  disorder  on  the  phonon  Raman  and  infrared  spectra.  Also  of  interest  is  the  time- 
dependence  of  phonon  and  electronic  Raman  spectra  of  quantum  wells  following  excitation  with  a 
picosecond  laser  pulse. 


SUMMARY  OF  RESEARCH: 

Our  Raman  study  of  phonons  in  GaAs/AsGaA*  super  lattices,  where  the  phonons  propagate  normal 
to  the  interfaces,  has  shown  that  the  longitudinal  acoustic  modes  propagate  through  the  interfaces  with 
only  small  zone-folding  effects.  On  the  other  hand,  the  optical  modes  show  strong  confinement  (much 
stronger  than  that  for  electrons  or  holes)  This  work  is  currently  being  extended  to  a  Raman  study  of 
.  phonons  propagating  parallel  to  the  interfaces. 

Raman  spectra  of  metastable  GaSb-Ge  alloys  grown  by  rf  sputtering  show  anomalously  broad  and 
asymmetric  peaks,  whereas  metastable  GaAs-GaSb  alloys  show  more  normal  behavior. 

Two  synchronously-pumped,  mode-locked  picosecond  dye  lasers  were  used  to  study  carriers  in 
undoped  GaAs-AlGaAs  quantum  wells  by  pulse-probe  Raman  scattering.  Preliminary  results  show 
that  the  higher  subband  levels  empty  in  about  50ps,  whereas  the  two  lowest  levels  take  several  hun¬ 
dred  picoseconds  to  empty. 

Multiple  and  single  quantum  wells  of  varying  well  sizes  were  investigated  in  detail.  Growth 
conditions,  source  type  was  optimized  which  led  to  the  achievement  of  extremely  sharp  transitions.  As 
a  result  identification  of  various  transitions  were  made  possible  for  the  first  time.  The  same  samples 
were  also  studied  under  hydrostatic  pressure,  under  high  excitation  pump  and  probe  configuration,  and 
in  the  presence  of  electric  field  both  parallel  and  perpendicular  to  the  interfaces.  In  addition  both  donor 
doped  and  acceptor  doped  samples  were  investigated.  The  quantum  wells  are  currently  being  investi¬ 
gated  under  extremely  low  cw  and  pulsed  excitation  conditions.  We  as  a  result  generated  a  wealth  of 
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information  regarding  not  only  the  origin  of  various  transitions  but  also  on  the  nature  of  crystal 
growth  and  island  formation. 
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SCIENTIFIC  OBJECTIVE: 

Metalorganic  chemical  vapor  deposition  (MOCVD)  is  a  sophisticated,  automated  epitaxial  growth 
process  for  thin,  high-quality  epitaxial  layers  of  GaAs  and  AlGaAs.  The  MOCVD  process  has  a  number 
of  advantages  important  for  the  development  of  complex,  thin-layer  quantum  limit  electronic  devices, 
such  as  high  uniformity  and  accurate  control  over  wide  ranges  of  thickness  (  <  20  A  to  10’s  of  /bim)  and 
composition  (  x  »  0  to  l)  and  interfaces  that  are  abrupt  on  a  monatomic  scale.  This  growth  process  has 
allowed  enormous  advances  in  the  development  of  heterostructure  lasers  and  other  optical  devices.  It  is 
the  objective  of  this  research  to  utilize  the  advantages  of  the  MOCVD  process  to  fabricate  high-quality 
epitaxial  electronic  devices  such  as  heterostructure  field  effect  transistors  (high  electron  mobility 
transistors,  HEMTs)  and  real-space  transferred  electron  devices. 


SUMMARY  OF  RESEARCH: 

The  research  has,  thus  far,  centered  primarily  on  forming  a  two-dimensional  electron  gas  at  the 
interface  of  an  undoped  GaAs-doped  AlGaAs  heterostructure  interface.  Considering  the  parts  of  this 
structure  individually,  we  have  found  that  there  is  a  large  degree  of  variability  in  available  starting 
source  materials,  such  as  arsine  gas  and  trimethylgallium.  These  sources  affect  the  background  doping 
level  and  maximum  electron  mobility  in  layers  of  undoped  GaAs.  Careful  choice  of  good  starting 
materials  and  optimization  of  growth  temperature  and  As-Ga  ratio  (which  affect  background  impurity 
incorporation)  has  allowed  growth  of  layers  with  background  mobilities  near  90,000  cm2/Vsec  (77  K). 
AlGaAs  doped  n-type  has  been  successfully  grown  over  the  entire  composition  range  with  controlled 
doping  and  reproducible  thickness.  Using  parameters  determined  from  a  model  previously  developed 
by  Hess  and  co-workers,  we  have  grown  single-interface  selectively-doped  heterostructures  that  show 
evidence  of  enhanced  mobility  in  the  conducting  notch  at  the  interface.  The  design  of  such  structures 
is  non-trivial  since  the  doping  and  thickness  relationship  must  be  such  that  the  wider-gap  layer  is  not 
completely  depleted,  yet  there  must  not  be  a  significant  parallel  (parasitic)  conduction  path.  This  can 
be  controlled  to  an  extent  by  using  gated  test  structures  and  optimizing  interface  conductance  by  vary¬ 
ing  gate  voltage. 


A  number  of  single-interface  structures  have  been  fabricated  into  large  area  gated  FETs.  The 
large  area  (1000  fim  source- to-drain,  970  fim  gate  width)  is  used  to  maintain  sufficiently  low  source- 
to-drain  electric  fields  for  Shubnikov-de  Haas  measurements.  Such  measurements  of  magneto-resistance 
as  a  function  of  field  strength  have  been  performed  in  a  super-conducting  magnet.  When  the  magnetic 
field  is  normal  to  the  two-dimensional  electron  gas,  oscillations  in  the  current-field  characteristics  are 
observed,  as  expected.  When  the  field  is  parallel  to  the  high  electron  density  interface,  only  bulk  mag¬ 
netoresistance  is  observed. 

In  summary,  enhanced  mobility  heterostructure  FETs  have  been  made  from  MOCVD  grown 
materials  which  have  reasonably  high  background  mobilities  and  exceptional  uniformity.  The  presence 
of  the  two-dimensional  electron  gas  has  been  confirmed  by  Shubnikov-de  Haas  measurements  on  large 
gated  structures. 


PUBLICATIONS: 


Journal  Articles: 

1.  T.C  Hsieh,  K.  Hess,  JJ.  Coleman,  and  PJD.  Dapkus,  "Carrier  density  distribution  in  modulation 
doped  GaAs  AlGaAs  quantum  well  heterostructures,"  Solid  State  Comm.  26,  1173  (1983). 
(JSEP,ARO) 

2.  D.  Widiger,  K.  Hess,  and  JJ.  Coleman,  "Two-dimensional  numerical  analysis  of  the  high  electron 
mobility  transistor,"  IEEE  Electron  Device  Lett.  EDL-5,  266  (1984).  (JSEP.ARO) 

3.  D.  Widiger,  LC  Kizilyalli,  K.  Hess,  and  JJ.  Coleman,  "Two-dimensional  transient  simulation  of 
the  high  electron  mobility  transistor,"  IEEE  Trans.  Elect.  Devices.  (JSEP.ARO) 


Con  ference  Pa  per  s : 


1.  D.  Widiger,  LC  Kizilyalli,  K.  Hess,  and  JJ.  Coleman,  "Two-dimensional  numerical  analysis  of  the 
high  electron  mobility  transistor,"  Proc.  First  Int.  Conf.  on  Superlattices,  Microstructures  and 
Microdevices  (Champaign,  1984),  p.  (JSEP.ARO) 

2.  D.  Widiger,  LC  Kizilyalli,  K.  Hess,  and  JJ.  Coleman,  "Two-dimensional  transient  simulation  of 
the  high  electron  mobility  transistor,"  Proc.  of  the  International  Electron  Devices  Meeting  (San 
Francisco,  1984),  p.  364.  (JSEPjkRO) 
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TITLE:  Studies  of  Transport  Phenomena  in  Semiconductors 


SENIOR  PRINCIPAL  INVESTIGATORS: 

K.  Hess,  Research  Professor 
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SCIENTIFIC  OBJECTIVE: 

This  research  involves  the  study  of  basic  properties  of  semiconductors,  semiconductor-heterolayers, 
new  device  concepts  and  device  modeling. 

Both  theoretical  and  experimental  methods  are  employed  in  each  of  these  categories.  We  are 
examining  a  variety  of  hot  electron  phenomena  and  their  effects  on  present  and  future  device 
performance,  especially  in  connection  with  modulation  doping.  The  experimental  studies  concern 
mainly  electronic  transport  in  heterolayers  in  high  electric  and  high  magnetic  fields. 

The  theoretical  studies  include  Monte  Carlo  simulations  of  electronic  transport  and  modeling  of 
ffl-V  compound  devices. 


SUMMARY  OF  RESEARCH: 

4.2.1.  Electronic  Transport  at  High  Energies  and  Transient  Phenomena* 

We  have  completed  a  study  of  high  field  transport  of  holes  in  GaAs  and  InP.  It  is  concluded  that 
steady  state  and  transient  transport  of  holes  in  these  materials  can  be  understood  in  significant  detail 
using  the  Monte  Carlo  simulation  technique.  We  have  been  able  to  predict  saturation  velocities  of  holes 
which  have  been  confirmed  by  recent  experiments  on  avalanche  photodiodes  at  Bell  Laboratories. 

Our  research  on  the  high  electron  mobility  transistor  (HEMT)  has  resulted  in  the  first  two- 
dimensional  device  model.  We  have  developed  both  a  Boltzmann  moment  equation  code  and  a  particle 
simulation.  The  model  has  proven  conclusively  the  importance  of  the  high  mobility  for  the  advantages 
due  to  a  low  source  access  resistance,  the  significance  of  velocity  overshoot  under  the  gate,  and  the 
importance  of  real  space  transfer. 

!  _ 

•  Thi*  work  w*i  supported  by  the  Joint  Services  Electronics  Program  under  contract  N00014-79-C-0424,  the  Artnv  Research 
Office  under  contract  DAAG-29-80-K-0069  and  the  Office  of  Naval  Research  under  contract  N00014-76-C-806. 
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A  third  project,  diffusion  in  very  high  electric  fields,  has  also  been  completed.  The  results  for  the 
anisotropic  diffusion  coefficient  are  important  for  the  assessment  of  the  significance  of  diffusion  in  high 
field  regions  of  devices  (MQS  transistors,  IMPATT  devices,  etc.). 

4.22.  Electronic  transport  in  semiconductor  heterolayer  structures. 

In  addition  to  work  on  the  HEMT,  reported  in  Section  4.2.1,  we  have  performed  research  on  two 
projects  related  to  heterolayer  structures.  A  model  has  been  developed  to  calculate  some  electronic  pro¬ 
perties  of  superlattices.  This  is  based  on  a  general  k.p  method  and  goes  beyond  the  naive  effective- 
mass-Kronig-Penney  approach.  It  has  been  applied  to  the  following  problems; 

(i)  We  have  studied  the  electronic  structure  of  superlattices  and  multiple  quantum  wells  under 
hydrostatic  pressure . 

(ii)  We  developed  a  realistic  model  to  study  the  optical  properties  and  the  dielectric  constant  of  super¬ 
lattices.  (For  details  see  publication  list.) 

A  second  project  has  been  an  in-depth  study  of  the  so-called  "Hickmott  experiment"  (J.  W.  Hick- 
mott  et  aL,  Phys.  Rev.  Lett.  52,  2053, 1984),  which  shows  that  oscillations  appear  in  the  I-V  characteris¬ 
tics  of  GaAs-AlGaAs  tunnel  junctions  in  the  presence  of  high  longitudinal  magnetic  fields.  Although 
magnetic  freeze-out  of  the  substrate  donors  is  necessary  for  the  onset  of  the  oscillations,  the  voltage 
period  is  rigorously  independent  of  the  magnetic  field  and  correspond  to  the  GaAs-L.0.  phonon  fre¬ 
quency.  We  have  shown  that  any  explanation  of  this  effect  by  ballistic  transport  has  to  be  incorrect 
under  the  given  experimental  conditions.  The  classic  explanation  on  phonon-assisted  tunneling  appears 
also  unsatisfactory  and  does  not  explain  several  experimental  observations.  We  propose  an  explanation 
based  on  space  charge  effects  induced  by  the  generation  of  local  perturbances  of  the  distribution  of  opti¬ 
cal  phonons.  This  interpretation  is  in  good  agreement  with  the  experimental  results  and  provides  a  gen¬ 
eral  basis  for  the  understanding  of  similar  effects  encountered  in  the  past  20  years  in  tunnel  junction 
structures. 

SIGNIFICANT  RESEARCH  ACCOMPLISHMENTS: 

The  most  significant  accomplishments  of  the  research  described  in  the  previous  sections  are,  in  our 
opinion,  the  explanation  of  the  Hickmott  experiment  and  the  development  of  the  Monte  Carlo  code  for 
the  high  electron  mobility  transistor.  This  code  represents  the  first  two-dimensional  model  of  its  kind 
and  is  currently  perceived  to  give  the  most  accurate  predictions  of  overshoot  effects. 

INTERACTIONS  AND/OR  TECHNOLOGY  TRANSFER: 

Monte  Carlo  codes  for  impact  ionization  phenomena  in  GaAs  and  Si  have  been  transferred  to 
Rockedyne,  Albuquerque,  in  connection  with  an  Air  Force  project  on  second  breakdown.  Some  of  the 
simulation  codes  for  transient  transport  have  been  developed  in  cooperation  with  the  US.  Army  Elec¬ 
tronics  Technology  and  Devices  Laboratory,  Ft.  Monmouth.  Interactions  on  heterolayer  transport  have 
continued  with  several  groups  in  the  Electronics  Technology  Division  of  the  Naval  Research  Labora¬ 
tory. 


PUBLICATIONS: 


Journal  Articles: 

1.  K.  Brennan  and  K.  Hess,  "Theory  of  high-field  transport  of  holes  in  GaAs  and  InP,"  Phys.  Rev.  B, 
VoL  29,  No.  10, 5581-5590  (1984).  (JSEP/ONR) 
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T.  Wang  and  K.  Hess,  "Calculation  of  high  field  diffusivity  by  a  many-  particle  Monte  Carlo 
simulation  including  a  complete  bandstructure  for  GaAs,"  /.  Appl.  Phys.,  October  (1984). 
(JSEP/ONR) 

P.  Gavrilovic,  J.  M.  Brown,  R.  W.  Kaliski,  N.  Holonyak,  K.  Hess,  M.  J.  Ludowise,  W.  T.  Dietz  and 
C  R.  Lewis,  "Resonant  tunneling  in  a  GaAs,  P  GaAs  strained-layer  quantum-well  heterost noc¬ 
ture,"  Solid  State  Comm.  52. 237-239  (3)  (1984).  1JSEP/ARO) 

T.  Wang  and  K.  Hess,  "Calculation  of  the  electron  velocity  distribution  in  high  electron  mobility 
transistors  using  an  ensemble  Monte  Carlo  method,"  J.  Appl.  Phys.  to  be  published. 
(JSEP/ONR/  ARO) 

J.  P.  Leburton,  "Size  effects  on  polar  optical  phonon  scattering  of  1-D  and  2-D  electron  gas  in  syn¬ 
thetic  semiconductors,"  J.  Appl.  Phys.  56,  2850  (1984).  (JSEP,  ARO) 

J.  P.  Leburton,  "Origin  of  the  current  oscillations  in  GaAs-AlGaAs  tunnel  junctions,"  Phys.  Rev. 
B  31  (1985).  (JSEP) 


PUBLICATIONS  (OTHER  SUPPORT): 

1.  D.  Widiger,  K.  Hess  and  J.  J.  Coleman,  "Two  dimensional  numerical  analysis  of  the  high  electron 
mobility  transistor,"  IEEE  Electron  Device  Lett*  VoL  EDL-5,  No.  7  (1984).  (ARO/ONR) 

2.  D.  Widiger,  L  C  Kizilyalli,  K.  Hess  and  J.  J.  Coleman,  "Two-dimensional  transient  simulation  of 
the  high  electron  mobility  transistor,"  Superlattices  and  Microstructures,  to  be  published. 
(ONR/ARO) 

3.  D.  Widiger,  L  C.  Kizilyalli,  K.  Hess  and  J.  J.  Coleman,  "Two-dimensional  numerical  analysis  of 
the  high  electron  mobility  transistor,"  IEEE  Transactions  on  Electron  Devices,  to  be  published. 
(ONR/ARO) 

4.  K.  Kahen,  J.  P.  Leburton  and  K.  Hess,  "General  model  of  the  transverse  dielectric  constant  of 
GaAs-ALAs  superlattices,"  Superlattices  and  Microstructures,  1  (1985).  (ARO/NASA) 

5.  K.  Brennan,  T.  Wang  and  K.  Hess,  "Theoretical  Calculation  of  Electron  Impact  Ionization  in 
Superlattice  Avalanche  Photodiodes,"  Applied  Phys.  Letters,  to  be  published.  (ONR/NSF) 

6.  T.  Wang,  K.  Hess  and  G.  J.  Iafrate,  "Time-dependent  ensemble  Monte  Carlo  simulation  for 
planar-doped  GaAs  structures,"  J.  Appl.  Phys.  to  be  published.  (ARO) 

7.  K.  Yokoyama  and  K.  Hess,  "Intersubband  phonon  overlap  integrals  for  AlGaAs/GaAs  single-well 
heterostructures.  Phys.  Rev.  B,  to  be  published.  (ARO) 

8.  J.  P.  Leburton  and  K.  E  Kahen,  "GaAs-AlGaAs  superlattice  band  structure  under  hydrostatic 
pressure:  an  analysis  based  on  the  envelope  function  approximation,"  Superlattices  and 
Microstructures  1, 49  (1985).  (ARO/NASA) 
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SCIENTIFIC  OBJECTIVES: 

The  primary  objective  of  this  research  program  is  to  investigate  energetic  particle-surface  interac¬ 
tions  which  control  the  nucleation  and  growth  kinetics,  chemistry,  and  physical  properties  of  alloy 
semiconductors  during  vapor  phase  crystal  growth  by  UHV  ion  beam  sputtering  and  accelerated-beam 
molecular  beam  epitaxy.  In  both  of  these  growth  techniques,  low  energy  ion-surface  interactions  allow 
an  efficient  coupling  of  kinetic  energy  to  the  growth  surface  upon  condensation  thereby  altering  the 
surface  reactivity  as  well  as  adsorption  and  adatom  diffusion  kinetics  allowing  single  crystal  film 
growth  at  lower  temperatures,  much  more  precise  control  over  dopant  incorporation  probabilities  and 
depth  distributions;  and  the  growth  of  unique  metastable  alloys.  This  work  is  being  pursued  from  both 
an  analytical  and  an  experimental  point  of  view  in  order  to  establish  a  detailed  understanding  of  fun¬ 
damental  film  growth  mechanisms. 

SUMMARY  OF  RESEARCH: 

Elemental  Incorporation  Probabilities  and  Depth  Distributions 

We  are  developing  a  general  model  for  predicting  elemental  incorporaton  probabilities  and  depth 
distributions  of  both  thermal  and  accelerated  species  incident  upon  the  growing  film  surface  during 
vapor  phase  deposition.  The  model  framework  contains  both  thermodynamic  and  kinetic  components 
and  accounts  for  dopant  diffusion,  surface  segregation,  and,  in  cases  where  the  dopant  flux  a  non- 
tbermal  (e.g,  "accelerated-beam"  MBE,  sputter  deposition,  and  plasma  assisted  CVDX  low  energy  ion 
bombardment  effects  including  implantation  and  enhanced  diffusion.  Initial  calculations  have  been  car¬ 
ried  out  for  GaAs  and  Si  MBE  in  which,  especially  for  Si,  common  dopants  are  typically  characterized 
by  strong  surface  segregation  and  temperature-dependent  incorporation  probabilities  cr. 
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The  model  predicts  that  a  critical  growth  temperature  exists  which  corresponds  to  the  transition 
from  "equilibrium"  segregrtion  to  kinetically-limited  dopant  segregation.  Moreover,  this  transition 
significantly  affects  both  tr  and  the  dopant  profile  broadening  as  a  function  of  the  film  growth  tempera¬ 
ture  T.  since  the  steady  state  segregation  ratio  r  (surface  to  bulk  dopant  fraction)  has  the  opposite 
dependence  on  T,  in  the  two  regimes.  In  addition,  we  have  predicted  for  the  first  time  that  because  of 
the  high  segregation  ratios  of  many  dopants  in  MBE  GaAs  and  Si  (r  can  easily  be  >  10s),  surface  struc¬ 
tural  phase  transitions  may  occur  during  film  growth  due  to  the  formation  of  ordered  dopant  over¬ 
layers.  This  is  expected  even  when  bulk  dopant  concentrations  are  <  101  cm-  .  We  have  recently 
obtained  the  first  experimental  evidence  for  this  effect  as  noted  below. 

Calculated  values  of  <ri T.)  and  calculated  depth  profiles  have  been  found  to  agree  very  well  with 
experimental  data  available  in  the  literature  for  both  group-m  acceptors  and  group-V  donors  in  Su 
"Anomalies"  in  previously  published  data  can  now  be  explained,  based  on  the  model,  as  being  due  to 
variations  in  dopant  desorption  rates  and/or  surface  segregation  kinetics  which  occur  in  response  to 
changes  in  film  growth  temperature  and  dopant  surface  coverage.  We  have  also  calculated  profiles 
resulting  from  arbitrarily  complex  doping  schedules  and  obtained  excellent  agreement  with  experimen¬ 
tal  results  for  Al,  Ga,  and  Sb.  In  addition,  the  model  predicts,  in  agreement  with  limited  experimental 
data,  that  a  growth  parameter  range  exists  in  which  abrupt  doping  profiles  can  be  obtained,  even  for 
dopants  which  exhibit  strong  surface  segregation.  We  have  recently  demonstrated  this  experimentally. 
Transition  temperatures  from  equilibrium  to  kinetically-limited  segregation  have  been  determined  for 
several  important  dopants  in  Si. 

We  have  carried  out  extensive  experiments,  in  collaboration  with  Prof.  J.-E.  Sundgren,  on  In  in 
Si.  In  is  a  deep  acceptor  and  from  model  calculations  is  expected  to  exhibit  strong  segregation  in  MBE 
Si.  The  In  incorporation  probability  crto  in  (100)  Si  grown  by  molecular  beam  epitaxy  was  found, 
using  secondary  ion  mass  spectrometry  (SIMS),  to  decrease  from  essentially  unity  at  film  growth  tem¬ 
peratures  T,  of  -  550*C  to  <  10-4  at  840  *C  SIMS  depth  profiles  of  both  uniformly-doped  and 
modulation-doped  samples  showed  evidence  of  strong  surface  segregation  with  the  amount  of  profile 
broadening  directly  related  to  <rta  (T,).  A  combination  of  in-situ  RHEED,  LEED,  and  Auger  electron 
spectroscopy  were  used  to  show  that  the  surface  segregation  rate  was  sufficient  over  a  wide  range  in  , 
and  In  to  Si  flux  ratios  to  cause  the  initial  (2xlXl00)  Si  surface  reconstruction  to  transform  to  (3x4; 
due  to  the  formation  of  an  ordered  In  surface  layer.  The  In  surface  coverage  in  the  (3x4)  state  was  ~ 
0.05-0.1  monolayer  even  though  the  bulk  In  concentration  was  <  2  x  10  cm  .  The  (2x1)  to  (3x4; 
surface  pha”»  transition  was  reversible  by  either  terminating  film  growth  and  re-evaporating  the  excess 
surface  In  or  terminating  the  In  flux  while  continuing  Si  film  growth. 

Initial  experiments  have  also  been  carried  out  to  determine  the  temperature  dependence  of  dopant 
depth  distributions  in  In-doped  (100)  Si  films  grown  by  MBE  The  films  were  departed  at  a  rate  of  ~ 
1  /«n  h-1  with  steady  state  In  concentrations  in  bulk  layers  ranging  from  1  x  10  to  6  x  10  cm  . 
The  amount  of  segregation-induced  profile  broadening  A|„  in  (2xlMl00)  films  was  found  from  SIMS 
analyses  to  reach  a  maximum  of  A,n  ~  600  nm  at  a  critical  temperature,  as  predicted,  and  to  decrease  at 
both  higher  and  lower  T, -values.  Calculated  In  depth  profiles  were  found  to  be  in  good  agreement  with 
experimental  results.  The  shape  of  calculated  profiles  was  determined  primarily  by  a  combination  of 
the  T.-dependent  In  incorporation  probability  and  the  steady  state  ratio  of  the  surface  to  bulk  In  frac¬ 
tions.  Indium  surface  coverages  ^  large  enough  to  be  observed  by  AES  were  obtained  in  layers  grown 
under  conditions  corresponding  to  strong  segregation. 

Experiments  are  presently  underway  to  measure  the  desorption  energy  Ej„  of  In  on  (100)  Si  as  a 
function  of  In  coverage  0ta,  to  investigate  In  overlayer  growth  mechanisms,  and  to  determine  the  sur¬ 
face  phase  diagram  for  In  on  Si  as  a  function  of  T,  and  0|„.  This  work  is  relevant  to  both  the  under¬ 
standing  of  dopant  segregation  during  MBE  Si  as  well  as  the  understanding  of  the  early  stages  ofthe 
growth  of  metal  films  on  Si.  A  variety  of  techniques  are  being  employed  including  AES,  LEED, 
RHEED,  and  modulate  beam  spectrometry  with  all  measurements  being  carried  out  in  UHV. 


.OV-Vv 

-V" 
.*•  A  A 


--'A. 

A'.'-'. 

■V-V-V--, 

*  V-  *.*. 


.V-V-'A 
■  •  •  «"  •  1 
»  U  *  ..  **  -  * 

av-v*;  v 

-  \  \  v 


13 


I  on/ Surface  Interactions 

The  choice  of  dopants  in  MBE  GaAs  and,  particularly,  in  MBE  Si  is  limited  by  the  fact  that  most 
common  dopants  (e.g.  Ga,  In,  As;  P,  and  Sb  in  Si  and  Zn,  Cd,  S,  Sc,  and  Te  in  GaAs)  have  very  low  ther¬ 
modynamic  sticking  probabilities.  However,  ionization  and  acceleration  of  these  dopants  to  relatively 
low  energies,  50  to  a  few  hundred  eV,  results  in  enhancements  in  net  incorporation  probabilities  <r  of 
up  to  8-10  orders  of  magnitude.  We  have  made  modifications  in  the  unique  single-grid,  solid-source, 
UHV  computable,  ion  gun  we  reported  last  year  (also  see  publication  list)  to  allow  us  to  form 
accelerated  molecular  beams  of  even  low  vapor  pressure  materials  such  as  In.  We  are  presently  investi¬ 
gating  the  incorporation  mechanism  of  accelerated  Zn  in  MBE  GaAs  and  accelerated  In  in  MBE  Si.  In 
both  cases  we  see  orders  of  magnitude  increases  in  <7  at  acceleration  energies  <  50  eV  corresponding  to 
extrapolated  ion  ’ranges"  of  the  order  of  or  less  than  a  monolayer.  Similar  experiments  carried  out 
using  accelerated  Sn,  where  <r  is  already  —  1  for  thermal  beams;  in  epitaxial  sputter-deposited  GaAs 
shows  evidence  of  a  large  decrease  in  the  dopant  segregation  rate  and  hence  a  much  more  uniform  dop¬ 
ing  profile. 

We  have  also  begun,  in  collaboration  with  Dr.  G.  Bajor  of  Harris  Semiconductor,  to  extend  our 
model  for  dopant/surface  (see  previous  section  and  publications  list)  interactions  to  include  ionized  and 
accelerated  beams.  We  have  succesfully  calculated  incorporation  probabilities  and  depth  profiles  of 
accelerated  ionized  dopants  in  MBE  GaAs  and  Si  and  used  the  model  to  explain  doping  behavior  in 
sputter-deposited  films.  In  the  latter  case,  preferential  sputtering  and  collisional  mixing  begin  to  play  a 
role  due  to  the  large  ion  fluxes  involved. 
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Growth  and  Characterization  of  Single  Crystal  Metastable  Semiconductors 

We  are  carrying  out  the  first  detailed  study  of  the  growth  and  physical  properties  of  new  single 
crystal  metastable  semiconductors.  The  key  feature  in  stabilizing  the  growth  of  these  materials  is  the 
controlled  use  of  low  energy  ion  bombardment  during  deposition  to  modify  elemental  sticking  proba¬ 
bilities  and  adatom  diffusivities  and  to  promote  dynamic  collisional  mixing  during  deposition.  We  have 
concentrated  primarily  on  the  study  of  (GaAs)1_xGex  because  of  the  importance  of  the  end-members; 
the  interest  in  Ge/GaAs  heterostructures;  and  the  fact  that  it  is  representative  of  a  new  subclass  of 
potentially  important  alloys;  (m— V^^IV^ 

Epitaxial  metastable  (GaAs)1_xGex  alloys  with  compositions  ranging  from  x  -  0  to  x  -  1  have  been 
grown  on  (100)  GaAs  substrates  by  ion  beam  sputtering  in  an  ultrahigh  vacuum  system.  Electron 
channeling,  double  crystal  X-ray  diffractometry,  and  X-ray  topography  analyses  indicate  that  the  films 
are  of  very  high  crystalline  perfection.  The  equlibirium  GaAs-Ge  pseudobinary  phase  diagram  has  been 
determined  by  differential  thermal  analysis  and  X-ray  diffraction/annealing  studies  to  be  a  simple 
eutectic  with  an  invariant  temperature  and  composition  of  880*C  and  ~  8  mole  %  GaAs.  We  have  also 
calculated  the  phase  diagram,  based  upon  a  modified  quasi-chemical  approach,  and  found  good  agree¬ 
ment  with  the  expeirmental  results.  The  model  is  now  being  extended  to  other  possible 
(m-V^/IV^  systems. 

Ether  n-type  or  p-type  conduction,  with  n  and  p  varying  over  several  orders  of  magnitude,  can 
be  obtained  by  varying  the  film  composition,  the  growth  temperature,  and  the  As  overpressure  during 
deposition.  The  As  incorporation  probability  <7*,  using  an  As4  beam  to  provide  the  overpresure;  was 
found  to  decrease  rapidly  with  increasing  Ge  concentration  in  the  films.  Detailed  measurements  of  o'Am 
using  both  Aa«  and  AS2  sources,  as  a  function  of  film  growth  temperature  and  Ge  concentration  are 
presently  being  carried  out  in  order  to  understand  and  model  surface  interactions  during  film  growth. 

Recent  high  resolution  triple-crystal  X-ray  diffractometry  measurements  of  (GaAs^.^GeA  on 
(100)  GaP  show  that  the  ratio  of  the  (200)  superlattice  to  (400)  fundmanetal  reflection  intensities 
decrease  precipitously  to  zero  near  x  *—  0.3.  This  is  the  first  direct  evidence  of  a  long-range 
order/disorder  phase  transition  in  these  materials.  It  is  in  agreement  with  our  earlier  prediction  based 
upon  a  model  developed,  in  collaboration  with  Professors  Jack  Dow  and  Kathie  Newman,  to  explain  the 
unusually  large  negative  bowing  in  the  direct  T-point  bandgap  E„  as  a  function  of  x.  We  are  presently 
carrying  out  XPS  measurements  to  map  out  the  valence  density  of  states  which  show  a  splitting  in  the  s 
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and  sp  peaks  due  to  the  transition. 

We  have  also  grown  epitaxial  metastable  alloys  of  (GaAsl^jCSnj),  and  Gej.gSn,  representing 
two  other  classes  of  metastable  semiconductors.  (GaAs^x^Sn^n  exhibits  essentially  no  solid  solubility 
and  requires  not  only  non-isovalent  substitution  on  both  cation  and  anion  sublattices  as  in 
(GaAs)(  1_x)(Ge2)»  but  non-isostructural  substitution  as  well.  That  is,  Sn,  at  temperatures  above  132*C  is 
a  tetragonal  semi-metal  and  the  growth  of  single  phase  (GaAs^x-x^Sn^,,  alloys  requires  the  stabiliza¬ 
tion  of  the  low  temperature  diamond  structure.  We  have  successfully  grown  alloys  which  are  not  only 
single  phase,  but  are  single  crystals,  on  (100)  GaAs  substrates,  with  Sn  concentrations  up  to  —  15%.  In 
addition,  high  temperature  phase  transition  studies  have  shown  that  under  the  proper  annealling  con¬ 
ditions,  (GaAs^j_x)(Sn2)x  can  be  decomposed  to  provide  dispersed  second-phase  regions  of  or- Sn  (dia¬ 
mond  structure)  which  are  completely  coherent  with  the  zinc-blende  lattice  at  temperatures  almost 
200*C  above  the  melting  point  of  (3— Sn!  Initial  Raman  results  show  that  the  Sn-based 
(m-V)(l_x)(IV2)x  alloys  exhibit  a  "two-mode"  behavior  while  high  resolution  X-ray  diffraction  analyses 
of  superlattice  and  fundamental  reflections  indicates  that  (GaAs^i_x)(Sn2)x  with  x  $  0.15  is  in  the 
ordered  zinc-blende  phase. 

Single  crystal  (GaAs)1_xGex  alloys  with  up  to  —  15  %  Sn  have  also  grown  on  (100)  GaAs  even 
though  the  equilibirium  Ge-Sn  phase  diagram  exhibits  essentially  no  mutual  solid  solubility  with  a 
eutectic  point  occurring  at  0.1  at  %  Ge  and  a  temperature  —  1*C  less  than  the  melting  point  of  pure  Sn, 
232*C.  This  alloy  represents  a  good  test  of  our  collisional  mixing  models  since  Sn  has  a  very  high 
diffusivity  in  Ge  and  hence  the  latitude  in  growth  parameter  space  (growth  temperature  T„  ion 
acceleration  potential  V„  and  deposition  rate  R)  is  extremely  narrow.  A  complete  growth  phase  map  as 
a  function  of  T„  V„  and  x  has  been  established.  For  a  given  value  of  x,  there  is  a  maximum  growth 
temperature  T,'  at  a  particular  value  of  V,  due  to  an  interplay  between  dynamic  collisional  mixing  of 
mobile  surface  adatoms  and  enhanced  diffusion  in  the  "bulk”  lattice.  The  critical  value  of  V,*  at  Tt* 
corresponds  to  an  ion  range  of  the  order  of  1  monolayer.  Both  T,‘  and  V,*  decrease  with  increasing  x. 

Phase  transition  studies  show  that  the  crystalline-metastable  to  equilibrium  transformation 
proceeds  through  a  much  different  reaction  path  than  for  (GaAs^1_x)(Gez)x.  There  are  no  intermediate 
metastable  states  in  the  (GaAs^^Ge,  due  to  structural  constraints.  The  earliest  observable  reaction  pro¬ 
duct  is  tetragonal  Sn.  Raman  results  indicate  a  one-mode  behavior  with  the  optical  phonon  frequency 
decreasing  rapidly  with  increasing  x. 

Optical  Information  Storage  Using  Explosive  Crystallization  of  In^GajSb  alloys 

This  work,  carried  out  in  collaboration  with  Dr.  CJE.  Wickersham,  who  received  his  PhD.  in  our 
group  under  JSEP  support  and  is  now  the  manager  of  the  Thin  Film  Group  at  Varian,  demonstrates  a 
practical  application  of  earlier  research  supported  by  JSEP  on  the  crystal  growth  of  ni-V  alloys  by 
sputter  deposition.  We  have  utilized  explosive  crystallization,  for  which  our  group  developed  the  first 
correct  model  in  1978,  of  amorphous  Inx^GaySb  for  optical  information  storage.  Information  encoding 
is  accomplished  using  a  low  power,  focused  laser  beam  to  crystallize  local  areas  of  the  film  by  normal 
nucleation  and  grain  growth.  The  remainder  of  the  film  is  "developed"  by  spontaneously  propagating 
an  explosive  crystallization  (EC)  front.  The  formation  can  be  retrieved  by  rastering  a  low  intensity 
laser  beam  across  the  storage  media  and  detecting  the  difference  in  optical  reflectivity  between  the 
roughened  EC  matrix  and  the  smooth  highly-reflecting  laser-annealed  regions.  A  15:1  signal-to-noise 
ratio  has  been  obtained  in  initial  experiments  with  5  pm  diameter  recorded  spots.  The  spot  temperature 
and  laser  power  densities  required  for  writing  with  this  approach  are  much  lower  than  in  comparable 
techniques  where  melting  or  vaporization  is  required. 

We  have  recently  established  experimentally  that  the  periodic  surface  undulations  in  films 
transformed  by  self -sustained  explosive  crystallization  are  due  to  thermal  instabilities  at  the  advancing 
amorphous/crystalline  boundary.  Transmission  and  scanning  electron  microscopy  studies  were  used  to 
show  that  the  amplitude  of  the  undulation,  which  can  be  as  large  as  micrometers,  is  directly  related  to 
the  transformation  temperature.  The  results  can  be  explained  using  instability  theory. 


SIGNIFICANT  RESEARCH  ACCOMPLISHMENTS: 

1.  Our  work  on  m-V  crystal  growth  by  sputter  deposition  has  resulted  in  the  development  of  an 
optical  storage  information  device  based  upon  explosive  crystallization  of  amorphous  Ini_,Ga^b. 
(See:  CJL  Wickersham,  G.  Bajor,  and  JE.  Greene,  Optical  Information  Storage  Using  Explosive 
Crystallization  in  Amorphous  Films,  J.  Vac.  ScL  TechnoL  Al,  18S7  (1983)  and  CE.  Wickersham, 
G.  Bajor,  and  JE.  Greene,  Temperature  Dependent  Formation  of  Surface  Undulations  in  Explo¬ 
sively  Crystallized  Films,  J.  Vac.  ScL  TechnoL,  in  press.) 

2.  Our  work  on  accelerated  beam  doping  during  MBE  has  resulted  in  the  design,  development,  and 
construction  of  a  unique  low-energy,  ultra-high  vacuum  compatible,  ion  source  which  can  be 
operated  with  gas;  liquid,  or  solid  source  material.  The  ion  gun  utilizes  single-grid  optics  and  has 
been  operated  with  Zn,  Ax  and  In  for  hundreds  of  hours  at  acceleration  emergies  as  low  as  20  eV 
(See:  A.  Rockett,  SA.  Barnett,  and  JE.  Greene,  A  Low-Energy,  Ultra-High  Vacuum,  Solid-Metal 
Ion  Source  for  Accelerated-Ion  Doping  During  Molecular  Beam  Epitaxy,  /.  Vac.  ScL  TechnoL  B2, 
306  (19841) 
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SCIENTIFIC  OBJECTIVE: 

Investigate  fundamental  plasma  and  chemistry  processes  in  a  type  of  magnetically  confined 
plasma  discharge  that  has  promise  for  reducing  radiation  damage  during  device  processing  by  dry  etch¬ 
ing. 


SUMMARY  OF  RESEARCH: 

This  project  is  designed  to  consist  of  essentially  three  one-year  phases.  The  first  phase  addresses  the 
basic  behavior  of  if  driven  magnetron  plasmas  and  the  resultant  plasma-surface  interactions.  It  calls 
for  the  assembly  of  an  apparatus  configured  to  permit  measurements  of  the  plasma  potential  distribu¬ 
tion  and  the  ion  bombardment  energy  at  surfaces  that  could  be  used  in  reactive  ion  etching  processes, 
and  for  the  making  of  such  measurements.  The  second  phase  would  explore  the  plasma  chemistry 
behavior  in  low  pressure  rf-driven  magnetron  discharges  using  CF^  and  CF4— 02  as  working  gases. 
The  third  phase  would  examine  silicon  etching  and  the  resultant  radiation  damage,  computer  modeling 
of  the  plasma  chemistry  process  and  the  extension  of  our  general  findings  to  plasma-assisted  deposition. 

Work  on  the  project  began  in  September  1984  with  the  addition  of  graduate  student  Guen  Yeom. 
The  primary  activities  during  this  first  six  month  period  have  been  in  preparing  a  laboratory  area  and 
in  designing  the  experimental  apparatus.  We  are  also  beginning  some  preliminary  computer  modeling 
of  the  chemistry  reactions  in  CF4H2  and  CF4— 02  plasmas  in  order  to  explore  the  systematic  changes 
which  can  occur  as  one  moves  from  the  gas  and  electron  density  conditions  that  are  typical  of  conven¬ 
tional  reactive  ion  etching  processes  to  those  conditions  of  low  gas  density  and  high  electron  density 
that  are  accessible  with  magnetrons.  A  third  activity  involves  plasma  modeling  studies  done  in  colla¬ 
boration  with  a  former  colleague  of  J.  Thornton’s.  The  objective  of  this  work  is  to  obtain  a  fundamen¬ 
tal  understanding  of  the  plasma  behavior  that  determines  the  current-voltage  relationships  in  glow 
discharges  and  therefore  that  influences  parameters  such  as  the  plasma  potential  which  are  important  to 
this  project. 

The  driving  force  in  our  work  is  the  rapid  growth  over  the  past  ten  years  or  so  in  the  extent  to 
which  rf  driven  glow  discharges  are  used  for  both  the  development  and  manufacture  of  new  materials 
and  devices.  Many  of  these  applications,  such  as  plasma-assisted  etching  and  plasma  assisted  chemical 
vapor  deposition,  involve  both  plasma-volume  and  plasma-surface  interactions.  Despite  the  fact  that 
glow  discharges  have  been  under  investigation  for  over  fifty  years,  even  the  general  form  for  the 
current-voltage  relationship  for  a  dc  planar  electrode  discharge  cannot  be  derived  from  first  principles. 
The  situation  is  even  more  complex  for  the  magnetron.  Often  the  effectiveness  of  a  particular  plasma- 
surface  interaction  depends  on  the  energy  of  positive  ions  incident  on  the  surface  in  question.  In  the  rf 
case  this  depends  on  the  time  variation  of  the  potential  of  the  surface  in  question  relative  to  the  plasma 
potential.  These  potential  variations  depend  in  turn  on  the  relative  areas  of  electrode  surfaces  that  are 


in  contact  with  the  rf  driven  plasma.  Only  very  recently  have  these  relationships  been  clarified  for  the 
relatively  simple  planar  diode  case.1  The  magnetron  case  is  considerably  more  complex  because  the 
magnetic  plasma  confinement  can  cause  the  effective  electrode  areas  to  be  much  different  from  the  phy¬ 
sical  areas.2  Furthermore,  and  perhaps  even  more  important,  the  relative  electron  and  ion  mobilities  in 
a  magnetically  confined  plasma  can  be  very  different  for  the  various  directions  that  connect  the  plasma 
to  the  electrodes. 

Therefore  our  Phase  One  experiments  are  designed  to  explore  the  relative  variations  in  plasma  and 
electrode  potential  that  develop  when  a  cylindrical  magnetron  is  driven  at  rf  frequencies.  The  funda¬ 
mental  glow  discharge  modeling  is  being  undertaken  in  an  attempt  to  pull  together  nearly  twenty 
years  of  collaborative  efforts  between  JA.  Thornton  and  A£.  Penfold3  in  order  to  develop  a  funda¬ 
mental  formulation  capable  of  explaining  the  current-voltage  behavior  of  both  conventional  and  mag¬ 
netron  discharges. 

The  experimental  apparatus  will  consist  of  an  20  cm  diameter  by  30  cm  high  cylindrical  vacuum 
chamber,  with  external  magnetic  field  coils  configured  to  provide  a  uniform  axial  magnetic  field.  The 
chamber  is  designed  to  mount  on  a  pumping  stand  with  an  oil  diffusion  pump,  and  to  accept  axially 
mounted  flange-type3  cylindrical  magnetron  electrodes  of  various  diameters.  Ports  are  available  for 
making  electrostatic  probe  measurements  in  both  the  radial  and  axial  directions  across  the  magnetron 
discharge,  and  for  making  optical  emission  spectroscopy  measurements.  The  initial  configuration  will 
mount  a  retarding  potential  electrostatic  analyzer4  on  the  grounded  chamber  wall.  It  is  anticipated  that 
as  the  Phase  One  experimental  work  proceeds,  an  electrostatic  analyzer  will  also  be  installed  within  one 
of  the  larger  cylindrical  magnetron  electrodes.  Power  supplies  are  available  for  driving  the  discharge 
at  frequencies  of  2  MHz  and  1X36  MHz.  Initial  experiments  will  consist  of  a  systematic  investigation 
of  the  variation  in  plasma  potential  (relative  to  the  electrode  surfaces)  as  a  function  of  the  relative  elec¬ 
trode  areas;  magnetic  field  strength,  and  rf  frequency,  while  using  inert  gases  such  as  He  and  Ar.  Par¬ 
ticular  attention  will  be  given  to  estimating  the  energy  of  the  ions  which  bombard  various  surfaces 
that  might  be  used  to  mount  wafers  for  plasma  etching. 

Our  preliminary  plasma  chemistry  calculations  treat  the  plasma  volume  as  a  thin  slab  with  uni¬ 
form  electron  density  and  temperature.  (It  is  assumed  that  the  electron  energy  distribution  function 
permits  the  use  of  an  effective  temperature.)  The  densities  of  various  active  species  are  calculated  by 
numerically  solving  a  reduced  set  of  rate  equations,  which  are  believed  to  account  for  the  major  reac¬ 
tions;  under  specified  conditions  of  working  gas  density  and  electron  density  and  temperature.  The  cal¬ 
culations  are  being  made  for  various  slab  thicknesses  and  plasma  conditions  that  span  from  those  used 
in  conventional  planar  diode  reactive  ion  etching  reactors  to  those  that  are  accessible  with  the  cylindri¬ 
cal  magnetrons.  Thus  gas  densities  are  being  varied  from  1013cm~3  to  lO^cm-3,  and  electron  densities 
from  109cm~3  to  10ucm~3,  with  an  effective  electron  temperature  of  3  eV.  The  calculations  are  being 
made  for  both  CF4— H2  and  CF4— 02  plasmas.  Our  objective  is  to  determine  the  effect  of  the  high  elec¬ 
tron  densities  and/or  low  gas  pressures,  which  are  available  in  magnetrons,  for  promoting  various  reac- 
-ion  paths  within  the  plasma  chemistry.  The  specific  goal  is  to  provide  guidance  for  planning  the  Phase 
Two  experiments.  The  calculations  will  be  refined  later  in  the  program  to  account  for  the  ionization 
balance  and  the  actual  nature  of  the  electron  distribution  function,  as  well  as  to  more  accurately  model 
the  geometry  of  the  cylindircal  magnetron  plasma. 
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SQENTIFIC  OBJECTIVE: 

6.0.  Introduction 

Electron  transport  is  accomplished  in  a  buried  channel  GaAs  structure  using  the  traveling  wave 
electric  field  of  a  Surface  Acoustic  Wave  (SAW)  generated  directly  in  the  GaAs.  The  electronic 
transport  over  long  distances  in  a  low  noise  environment  makes  this  technology  highly  dependent  on 
GaAs  layer  parameters. 


6.1.  Charge  in  Gallium  Arsenide  Layers  Transported  by  Traveling  Wave  Electric  Fields 

A  new  type  of  charge  transfer  based  on  a  non-linear  interaction  between  an  electric  field 
traveling  waves  and  electrons  in  Gallium  Arsenide  has  recently  been  reported.1  In  this  structure, 
buried  channel  electron  transport  is  accomplished  by  using  a  traveling  wave  electric  field  generated 
directly  in  Gallium  Arsenide.  A  theoretical  investigation  of  material  on  the  basis  of  transport 
efficiencies  is  possible,  primarily  because  the  traveling  wave  excitation  provides  a  nearly  uniform  drift 
field  so  that  carrier  transport  occurs  entirely  under  the  influence  of  continuous  electric  fields. 

Significant  progress  has  been  made  towards  the  experimental  demonstration  of  the  theoretical 
transport  efficiency.  The  charge  transport  efficiency  is  believed  to  be  primarily  material  limited.  A 
new  simplified  epitaxial  layer  structure  is  used  to  form  the  buried  channel,  and  it  is  grown  directly  on 
a  GaAs  semi-insulating  substrate.  This  simplified  layer  structure  shows  a  lower  defect  density  and  a 
higher  transport  efficiency.  In  addition,  the  mesa  structures  previously  used  for  lateral  confinement  of 
the  charge  transport  have  been  replaced  by  a  Schottky  guard  rail  structure.  These  configuration 
improvements  resulted  in  a  transfer  efficiency  in  excess  of  .996  when  the  traveling  wave  frequency  was 
358  MHz. 

The  transfer  efficiency  due  to  proportional  loss  was  determined  by  leading  edge  measurements  and 
measured  as  a  function  of  signal  leveL  These  data  compared  closely  to  the  theory  predicted  diffusion- 
limited  efficiency  for  traveling  wave  electric  field  strengths  of  8,000  volts/centimeter. 


At  large  signal  strengths,  the  ballooning  effect  observed  due  to  the  diffusion  limited  transport 
followed  the  predictions  of  the  theoretical  model.  However,  the  exact  cause  of  the  transfer  inefficiency 
at  smaller  charge  loads  is  not  precisely  known.  This  inefficiency  is  believed  to  be  due  to  static  channel 
potential  variations  arising  from  structures  at  the  surface  and  to  trapping  sites  in  the  Gallium  Arsenide 
layers. 

The  depth  of  buried  channel  transport  is  adjusted  by  applying  bias  to  a  Schottky  plate  located  on 
top  of  the  buried  channel.  It  has  been  clearly  demonstrated  that  the  small  charge  load  transport 
inefficiency  is  directly  related  to  the  position  of  the  charge  transport  in  the  channel,  possibly  due  to  a 
variation  in  trapping  site  density  in  terms  of  depth. 

The  transport  efficiency  of  electrons  transported  by  traveling  wave  electric  fields  appears  to  be  a 
very  sensitive  means  for  measuring  trapping  densities  and  lifetimes.  Further  work  should  be  pursued 
to  determine  the  value  of  these  data  in  the  characterization  of  Gallium  Arsenide  epitaxial  layers. 
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SCIENTIFIC  OBJECTIVE: 

The  goal  of  this  research  is  a  better  understanding  of  the  InGaAsP/InP  quaternary  alloy  system 
for  application  to  heterostructures  and  heterostructures  devices,  particularly  avalanche  photodiodes  and 
microwave  devices.  In  many  of  these  applications  it  is  important  to  have  the  capability  of  growing 
high  purity  material  with  abrupt  compositional  changes  at  heterojunction  boundaries.  Because  of  this, 
the  purity  attainable  with  the  hydride  growth  technique  will  be  studied  and  compared  with  the 
results  obtained  from  other  epitaxial  growth  techniques  such  as  MOCVD  and  MBE  In  this  work, 
undoped  epitaxial  layers  of  InP,  GaAs,  and  InGaAs  are  grown  using  the  hydride  technique,  and  the 
resulting  samples  as  well  as  those  obtained  from  other  laboratories  are  characterized  using  as  a  wide 
variety  of  techniques.  These  techniques  include  Hall  coefficient  measurements  and  analyses;  low  tem¬ 
perature  photoluminescence  measurements,  far  infrared  photothermal  ionization  measurements,  double 
crystal  x-ray  diffraction  I-V  and  microwave  time-of -flight  techniques.  In  addition,  the  detector  struc¬ 
tures  fabricated  in  this  work  are  characterized  using  relative  spectral  response  measurements,  I-V  meas¬ 
urements,  photomultiplication  measurements  and  related  techniques  to  help  characterize  the  InGaAsP 
materials  used  in  these  structures.  Correlation  of  the  growth  parameters  with  the  results  of  these  meas¬ 
urements  will  hopefully  lead  to  optimization  and  better  control  of  the  epitaxial  growth  process  and  an 
understanding  of  the  physics  of  transport  processes  in  these  materials  which  is  important  for  develop¬ 
ment  of  higher  performance  devices. 


SUMMARY  OF  RESEARCH: 

High  purity  epitaxial  layers  are  critical  for  many  types  of  high  performance  m-V  devices.  Buffer 
layers  in  microwave  FETs,  multiplication  regions  in  APDs,  and  the  drift  regions  of  IMPATT  all  require 
low  net  carrier  concentrations  for  optimum  device  performance.  With  the  inevitable  small  fluctuations 
in  Nd  and  NA  the  only  way  to  reproducibly  maintain  a  low  net  concentration,  ND— NA,  is  to  grow  high 
purity  material.  When  applied  with  proper  caution,  the  77  K  Hall  mobility,  /sm^K,  provides  a  good 
indication  of  sample  purity,  and  is  the  usual  figure  of  merit  for  high  purity  IB-V  material. 


The  Hall  mobilities  and  carrier  concentrations  for  eight  of  the  best  GaAs  epitaxial  layers  grown  in 
the  hydride  VPE  system  are  presented  in  Table  2.1.  quality  and  reproducibility  obtained.  These  layers 
were  grown  all  at  a  GaCl  partial  pressure,  P<j«a.  of  10-2  atom  with  the  noted  D2-V  ratios.  The  parame¬ 
ter  which  we  have  found  to  have  the  greatest  effect  on  the  impurity  incorporation  is  the  ID-V  ratio  for 
the  range  of  growth  conditions  attempted.  In  these  experiments  the  m-V  ratio  was  varied  by  adjusting 
the  arsine  flow  while  maintaining  P^q,  the  growth  and  source  zone  temperatures  and  total  flow  in  the 
reactor.  The  effect  of  this  variation  in  PAlHj  is  illustrated  in  Fig.  2.2.  There  are  two  opposing  effects  to 
be  noted  in  this  set  of  photothermal  ionization  spectra.  The  middle  peak  due  to  sulfur  is  nearly  con¬ 
stant  for  the  range  considered  while  Si  concentration  decreases  and  Ge  concentration  increases  as  AsH3 
increases.  Both  Si  donors  and  acceptors  (the  acceptors  are  determined  from  PL)  are  suppressed  by  adding 
AsH>  The  Si  is  predominantly  due  to  SiCl*  which  due  to  the  chlorine  atoms  couples  it  to  the  growth 
reaction.  Because  both  reactions  liberate  HC1  the  additional  increase  in  AsH3  suppresses  Si.  Germanium 
donors  and  acceptors  increase  with  AsH«  flow.  Arsine  flow  was  the  only  growth  parameter  found  to 
influence  Ge  contamination.  This  along  with  the  lack  of  Ge  in  InP  grown  in  the  same  reactor  indicates 
that  the  arsine  tank  itself  is  the  Ge  source  in  this  system.  IILV  ratio  studies  in  MOCVD  reactors  also 
show  similar  evidence  for  Ge  in  semiconductor  grade  arsine.  If  the  Ge  were  a  hydride  (germane)  it 
would  not  be  suppressed  by  the  growth  reaction’s  liberation  of  HCL  Because  it  has  the  same  molecular 
mass  as  AsH«  it  is  difficult  to  detect  by  the  means  used  to  analyze  gases  in  the  part  per  billion  range  this 
mid  1014  cm-3  doping  level  would  indicate. 

Although  the  addition  of  oxygen  in  both  hydride  and  chloride  reactors  had  been  shown  to 
decrease  the  concentration  of  donors  and  acceptors,  the  evidence  was  from  Hall  effect  measurements  at 
room  and  liquid  nitrogen  temperatures.  This  led  to  some  ambiguity  about  whether  the  phenomenon 
was  a  decrease  in  the  activity  of  Si,  a^,  or  a  result  of  the  oxygen  being  incorporated  as  a  deep  donor  or 
trap.  The  photothermal  ionization  spectra  for  three  samples  grown  with  the  indicated  added  oxygen 
and  otherwise  identical  growth  conditions  in  Fig.  1  illustrate  that  it  is  indeed  the  decrease  in  ay  which 
is  responsible  for  the  decrease  in  Nq.  This  study  has  also  provided  some  additional  insight  into  the 
kinetics  of  the  model  used  to  explain  the  reduction  of  Si.  In  this  model  the  H2  and  HC1  reduce  the  hot 
silica  glass  walls  of  the  reactor  forming  volatile  chlorosilanes  which  are  transported  to  the  growth  zone 
and  incorporated  into  the  growing  crystal.  From  thermodynamic  arguments  alone  it  is  possible  to  show 
the  decrease  in  a^  and  02  but  these  considerations  do  not  account  for  the  absence  of  the  suppression  of 
when  the  02  is  added  downstream  of  the  gallium  source.  If  the  component  gas  species  came  to  equili¬ 
brium  at  each  step  as  the  usual  thermodynamic  arguments  assume,  there  would  be  no  difference.  That 
the  suppression  only  occurs  when  02  is  added  to  the  source  zone  indicates  the  formation  of  a  chlarosi- 
lane  is  the  step  in  the  process  affected  by  the  addition  of  0?.  Once  they  are  formed  they  are  too  stable  to 
equilibrate  with  O2  introduced  afterward  in  the  short  time  they  remain  in  the  reactor. 
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WORK  UNIT  NUMBER  8 


TITLE:  Direct  Examination  of  the  Metal-Semiconductor  Interface 


SENIOR  PRINCIPAL  INVESTIGATOR: 


Gert  Ehrlich,  Research  Professor 


SCIENTIFIC  PERSONNEL  AND  TITLES: 


J.  D.  Wrigley,  Research  Associate 


SCIENTIFIC  OBJECTIVE: 


Silicon  interfaces  have,  for  a  long  time,  been  of  paramount  importance  in  device  technology. 
However,  the  atomic  events  significant  for  the  formation  of  metal-silicon  interfaces  or  for  their  stabil¬ 
ity  are  not  well  understood.  The  aim  of  this  project  is  twofold:  (l)  to  examine  the  interactions  of  sili¬ 
con  atoms  with  metal  surfaces,  using  field  ion  microscopic  techniques  to  attain  a  direct  view,  on  the 
atomic  level,  of  the  processes  important  in  layer  formation;  and  (2)  to  characterize  the  atomic  structure 
of  the  interface.  Particular  emphasis  will  be  upon  the  variations  in  behavior  dictated  by  the  structure 
of  the  substrate,  as  well  as  upon  the  effect  of  foreign  atoms  upon  the  properties  of  the  layers  formed. 


SUMMARY  OF  RESEARCH: 


We  have  concentrated  our  efforts  on  characterizing  interactions  between  a  single  silicon  atom  and 
a  metal  atom  at  a  surface,  so  as  to  provide  information  about  the  energetics  of  silicon-metal  overlayers. 
Recent  work  has  emphasized  tungsten  and  rhenium  complexes  with  silicon  on  the  (110)  plane  of 
tungsten.  In  order  to  establish  the  strength  of  interactions  between  silicon  and  metal,  quantitative  stu¬ 
dies  of  the  kinetics  of  dissociation  of  W-Si  have  been  made  at  different  temperatures;  and  are  shown  in 
Fig.  1.  These  led  to  a  dissociation  energy  of  0.85eV,  which  is  just  the  same  as  the  barrier  to  motion  of  a 
single  silicon  on  W(110X  We  conclude  that  the  apparent  binding  energy  between  silicon  and  tungsten 
is  small,  and  within  the  limit  of  error  of  the  measurements.  This  is  quite  different  from  Pd-Si  on 
W(llO),  where  our  previous  studies  indicated  a  binding  energy  of  at  most  0.26  eV  and  suggests  very 
strong  chemical  effects. 

In  the  last  few  months  we  have  made  an  important  discovery:  it  is  possible  to  distinguish  the 
location  of  the  silicon  from  that  of  the  metal  atom  in  a  cluster.  With  this  entirely  new  capability  we 
have  been  able  to  examine  the  mechanism  of  motion  of  metal-silicon  complexes.  For  W-Si,  movement 
appears  to  occur  by  the  jump  of  a  single  atom  to  form  an  intermediate,  shown  in  Fig.  2,  in  which  both 
atoms  are  aligned  along  [001  ].  An  alternative  intermediate,  oriented  along  [llOl  has  a  larger 
interatomic  spacing,  and  therefore  presumably  a  higher  energy.  No  indication  has  been  found  that  this 
contributes  to  the  motion  of  W-Si  below  280K.  For  Re-Si,  however,  both  diffusion  paths  are  possible, 
even  though  transitions  through  the  [001]  intermediate  dominate.  Quantitative  studies  of  the  different 
diffusion  processes  of  Re-Si  on  W(110)  have  been  carried  out,  allowing  us  to  derive  the  rates  of  the 
different  types  of  atomic  jumps,  as  well  as  the  barriers  over  which  they  occur.  We  are  now  in  the 
process  of  deriving  information  about  the  silicon-metal  potentials  from  these  measurements. 
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Hard  sphere  models  of  hetero-dimer  diffusion  on  W(110).  At  left,  movement  of  grey  ada¬ 
tom  in  the  dimer  occurs  via  a  [001]  oriented  intermediate.  Even  after  diffusion,  grey  atom  is 
always  at  right  of  black  atom.  In  sequence  at  right  a  jump  of  grey  adatom  leads  to  a  [lTO] 
intermediate,  and  finally  to  a  dimer  with  adatoms  inverted  compared  with  final  state  at  left. 
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Conference  Papers: 

1.  John  D.  Wrigley,  Best  Poster  Award,  American  Vacuum  Society  Meeting,  Reno,  Nevada,  Dec. 
1984.  (JSEP) 

ADDITIONAL  INFORMATION: 

Awards  and  Honors 
Gert  Ehrlich 


Guggenheim  Fellow,  1984-85 


WORK  UNIT  NUMBER  10 
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SCIENTIFIC  OBJECTIVE: 

The  goal  of  this  research  is  to  understand  the  basic  mechanism  of  charge-density  wave  (CDW) 
transport  in  the  "one-dimensional"  metals  NbSe>  TaS>  and  related  materials.  This  CDW  conduction 
represents  the  second  known  example  of  current  flow  due  to  motion  of  a  macroscopically  occupied 
quantum  state,  superconductivity  being  the  first.  John  Bardeen  proposed  that  CDW  acceleration  in  an 
electric  field  takes  place  via  a  collective  tunneling  process  [l]  and  that  the  ac  response  of  these  systems 
can  be  characterized  by  applying  the  photon-assisted  tunneling  theory  originally  developed  by  the 
principal  investigator  for  superconductor-insulator-superconductor  (SIS)  millimeter  wave  mixers:  [2,3] 
This  controversial  hypothesis  implies  that  quantum  phenomena  are  observed  in  CDW  systems  in  the 
MHz  frequency  region,  where  hu>  *s  10-6  kT  at  temperatures  as  high  as  T  -  200  K-  Such  behavior  is 
theoretically  possible  because  very  large  numbers  of  electrons  are  acting  together,  so  that  the  total 
energy  involved  is  larger  than  thermal  energies.  If  this  interpretation  is  correct,  CDW  materials  could 
become  important  systems  for  both  long-wavelength  detectors  and  for  basic  studies  of  macroscopic 
quantum  tunneling  and  coherence.  An  extensive  experimental  program  is  now  underway  in  our 
laboratory  to  characterize  the  dynamics  of  CDW  response,  and  to  test  hypotheses  based  on  the  tunneling 
model. 


SUMMARY  OF  RESEARCH: 

During  the  past  year,  experimental  work  has  focused  on  developing  a  phenomenological  model 
which  describes  the  complete  response  of  NbSe>  TaS*  and  related  CDW  conductors  to  arbitrary  time- 
dependent  applied  potentials.  Our  previous  studies  [4,5]  had  demonstrated  that  the  nonlinear  conduc¬ 
tion  and  mixing  induced  by  small  ac  signals  could  be  successfully  interpreted  using  the  tunneling 
theory  when  the  applied  dc  voltage  was  above  the  threshold  value.  Recent  measurements  of  the  non¬ 
linear  response  below  threshold,  together  with  experiments  involving  large  as  well  as  small  ac  signals, 
have  led  to  a  relatively  simple  model  capable  of  characterizing  the  full  range  of  experimental  observa¬ 
tions  in  the  MHz  frequency  region.  [6,7]  The  model  is  based  on  the  tunneling  theory,  and  we  believe 
that  this  consistent  and  quantitative  interpretation  of  the  complex  field  and  frequency  dependent  con¬ 
duction  in  these  systems  is  powerful  evidence  in  favor  of  the  quantum  tunneling  theory  of  CDW 


I 


SIGNIFICANT  RESEARCH  ACCOMPLISHMENTS: 

During  the  past  two  years,  work  on  the  new  ultra-low  noise  millimeter  wave  receivers  which 
utilize  photon-assisted  tunneling  in  superconductor-insulator-superconductor  (SIS)  junctions  has  pro¬ 
gressed  rapidly.  Many  new  programs  have  been  established  to  engineer  SIS  receivers  in  the  region 
between  36  GHz  and  240  GHz  in  laboratories  around  the  world,  these  nonlinear  tunneling  devices  are 
theoretically  capable  of  quantum  limited  receiver  temperatures  TjpL  =  hv/k,  and  a  practical  receiver 
system  was  recently  constructed  [8]  with  a  single-sideband  noise  temperature  TR  -  68*C  K  *=  10  h*»/k 
operating  at  115  GHz.  Much  lower  receiver  temperatures  should  be  realized  during  the  next  several 
years,  and  this  revolutionary  advance  in  detection  capability  is  expected  to  have  a  major  impact  on 
radio  astronomy  research.  Future  intelligence  and  communication  applications  in  space  will  also’ be  pos¬ 
sible.  A  complete  review  of  the  status  of  these  receivers  entitled  "Quantum  detection  at  millimeter 
wavelengths"  has  been  written  in  collaboration  with  MJ.  Feldman,  and  will  be  published  in  Reviews 
of  Modem  Physics. 
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SCIENTIFIC  OBJECTIVE: 

The  general  objective  of  this  research  unit  is  to  study  the  interaction  of  excited  molecules,  atoms 
or  molecular  fragments  with  solid  surfaces  with  particular  emphasis  on  the  etching  or  growth  of 
semiconducting  materials.  The  excited  species  may  be  produced  by  a  discharge,  electron  beam, 
ultraviolet  or  visible  photons;  or  combinations  of  the  three.  Inasmuch  as  these  excitation  means  produce 
an  environment  which  is  far  from  thermodynamic  equilibrium,  one  can  expect  processes  which  deviate 
considerably  from  that  found  in  conventional  chemical  processing.  Our  immediate  experimental  goals 
are  to  utilize  many  of  the  modern  sophisticated  diagnostic  tools;  such  as  laser  induced  fluorescence, 
optogalvanic  spectroscopy,  mass  spectroscopy,  and  Raman  spectroscopy  to  elucidate  the  physical 
processes  occurring  in  the  volume  adjacent  to  the  surface  being  etched  or  grown.  One  further  goal  is  to 
grow  (by  laser  or  discharge  techniques)  electronic  materials  and  structures  which  are  difficult  to 
fabricate  by  conventional  means. 

SUMMARY  OF  RESEARCH: 

A.  Laser  Growth  of  Films 

During  this  past  year,  emphasis  has  been  placed  on  examining  the  effect  of  illuminating  surfaces 
with  visible  or  UV  laser  radiation  during  the  CVD  growth  of  films.  While  this  work  is  still  in  an  ini¬ 
tial  phase,  it  has  been  observed  that  even  low  fluences  of  laser  radiation  significantly  alter  the  film 
growth  rate  and  quality.  For  example,  for  Ge  films  grown  by  the  pyrolysis  (CVD)  of  GeEL  at  tempera¬ 
tures  just  beyond  300*C,  the  growth  rate  increases  by  roughly  an  order  of  magnitude  when  the  sub¬ 
strate  is  illuminated  with  ArF  laser  (X  -  193  nm)  radiation.  Preliminary  measurements  indicate  a 
dramatic  change  in  the  index  of  refraction  of  the  deposited  film  as  well.  Similar  observations  have 
been  made  when  an  Ar*  ion  laser  beam  (X  -  514  nm)  strikes  the  surface.  In  the  next  year,  we  expect  to 
explore  the  laser  wavelength  and  intensity  dependence  of  this  effect  as  it  shows  considerable  promise 
for  permitting  the  growth  of  high  quality  films  at  significantly  reduced  substrate  temperatures. 

We  have  also  recently  demonstrated  the  feasibility  of  detecting  minute  impurities  in  gases  by 
laser-induced  breakdown  (LIB).  A  pulsed  visible  laser  is  focused  into  a  cell  containing  a  small  sample 
of  the  gas  in  question  and  emission  resulting  from  UB  in  the  gas  is  dispersed  (spectrally)  and  detected. 
Already,  argon  has  been  detected  in  GeH«  at  the  sub- ppm  level  without  the  need  for  synchronous 
detection.  Therefore,  this  technique  shows  promise  for  detecting  impurities  which  cannot  be  studied  by 


conventional  means  such  as  mass  spectroscopy. 

.Finally,  in  the  past  year,  the  spectroscopic  studies  of  a  reactor  in  which  germanium  films  are 
grown  by  photolytic,  laser-induced  CVD  were  completed  and  the  results  will  be  published  in  the  May, 
1985  issue  of  the  Journal  of  Applied  Physics.  The  results  clearly  point  to  the  roles  played  by  GeH2 
and  GeH3  in  the  deposition  process. 


B.  Discharge  Etching 

Ip  the  last  progress  report,  it  was  mentioned  that  considerable  insight  could  be  obtained  on  the 
kinetics  involved  in  discharges  used  in  plasma  processing  by  pulsing  the  discharges  and  observing  the 
time  evolution  of  the  plasma  parameters  such  as:  optical  emission,  dissociation,  and  the  formation  of 
new  species.  The  kinetics  of  the  NF3  system  has  been  examined  in  greatest  detail  and  illustrates  the 
wealth  of  information  available  from  the  pulsed  discharge. 

A  discharge  can  be  initiated  on  a  time  scale  of  a  few  microseconds  or  so,  whereas  the  approach  to 
equilibrium  of  any  of  those  parameters  requires  hundreds  of  microseconds.  For  instance,  consider  the 
"well  known"  optical  emission  from  fluorine  at  704  nm  from  etching  discharges.  In  our  pulse  discharge, 
the  optical  intensity  would  first  grow  linearly  with  time  and  then  saturate  at  its  equilibrium  value. 

All  of  our  experiments  are  consistent  with  a  two-step  process;  the  first  electron  to  dissociate  the 
NF3  and  the  second  to  excite  it  to  the  state  from  which  the  radiation  is  emitted. 


e  +  NF3  -*  NF^  +  xF  +  e 

(1) 

e  +  F-F  +  3 

(2) 

F  -  F  +  h  v  (704nm) 

(3) 

Note  that  this  scenario  is  in  contrast  to  that  proposed  by  Donnelly,  Flamn,  et  ai  [l]  who  suggested 
that  dissociative  attachment  might  be  a  major  source  of  atomic  fluorine.  However,  the  square  law 
dependence  on  electron  density  (or  current)  implied  by  Eqns.  1  and  2  was  easily  documented,  and  thus 
a  dissociative  attachment  sequence  Eq.  4-6  would  require  a  simultaneous  detaching  and  excitation  colli¬ 
sion 

e  +  NF3  -*  NF2  +  F~  (4) 

e  +  F"  -  F*  +  2e  (5) 

F*  -*  F  +  h?  (704nm)  (6) 

Previously,  [2]  we  had  demonstrated  that  100  mJ  of  excimer  laser  radiation  at  308  nm  was 
sufficient  to  completely  deplete  the  negative  ion  population  (F“)  and  thus  interrupt  this  kinetic 
sequence.  No  such  interruption  was  observed  and  thus  the  channel  indicated  by  Eqns.  4-6  was  not  con¬ 
sidered  important. 

The  channel  leading  to  the  excitation  of  the  C-B  band  of  N>  a  specie  which  is  formed  in  the 
discharge,  was  also  identified  and  is  represented  by  the  following  sequence. 

e  +  NF3  -  NF  +  2F  +  e  (7) 

(8) 
(9) 
(10) 


NF  +  NF  -»  N2  +  2F 
e  +  N2  -*  N2  +  e 
Nj  N2  -  +h  (337nm) 


It  takes  one  electron  to  produce  each  NF  molecule  and  thus  the  ground  state  population  on  N2 
would  vary  as  the  square  of  the  electron  density.  The  final  electron  would  excite  the  molecule  to  the  C 
state  from  which  radiation  was  observed.  This  sequence  was  verified,  and  thus  scenario  indicated  by  7 
and  8  appears  to  be  the  logical  path  for  the  formation  of 

It  appears  that  the  approach  to  equilibrium  of,  say,  the  F  emission  is  primarily  due  to  the  deple¬ 
tion  of  the  NF3  donor  gas.  Thus  by  utilizing  the  time  history  of  the  optical  emission,  one  can  determine 
the  appropriate  rate  coefficients  to  predict  the  percentage  dissociation  in  the  discharge.  Excellent  agree¬ 
ment  between  the  predictions  and  the  experiments  has  been  achieved.  This  work  has  been  described  in 
greater  detail  in  Ref.  3. 

In  the  course  of  this  work,  it  was  observed  that  the  addition  of  N2  to  the  NF3  discharge  increased 
the  etch  rate  by  as  much  as  509b.  Although  the  effect  was  established,  the  mechanism(s)  for  this 
enhancement  has  not  been  identified.  Some  of  the  possibilities  under  consideration  are:  (a)  a  chemical 
reaction  with  some  of  the  free  radicals  NFy_x  liberating  excess  fluorine:  (b)  An  electronic  energy 
transfer  from  the  nitrogen  metastable  N2  (A)  to  the  donor  gas;  or  (c)  a  change  in  the  electron  distribu¬ 
tion  function. 


INTERACTIONS  AND/OR  TECHNOLOGY  TRANSFER: 

The  research  described  above  and  in  Ref.  3  has  been  transferred  en  masse  to  an  Air  Force  contract 
involving  the  deposition  of  silicon  from  helium  -  SiH«  mixtures. 
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SCIENTIFIC  OBJECTIVE: 

There  are  four  principal  objectives  of  this  effort.  The  first  is  to  develop  efficient  iterative  tech¬ 
niques  for  solving  the  problems  of  electromagnetic  scattering  from  the  coupling  into  complex  struc¬ 
tures:  The  second  is  to  study  techniques  for  reducing  the  radar  cross-section  of  targets.  The  third  is  to 
analyze  frequency  selective  surfaces  (FSS)  for  radomes  and  other  antenna  applications.  Finally,  the 
fourth  objective  is  to  investigate  the  interconnection  problem  in  Packaging  and  develop  CAD  tools  for 
modeling  the  cross-talk,  coupling,  and  pulse  distortion  in  high  speed  digital  circuits. 


SUMMARY  OF  RESEARCH: 

In  this  project  we  have  been  concentrating  on  four  different  aspects  of  electromagnetic  scattering 
and  coupling  problems.  The  first  of  these  is  aimed  toward  developing  analytical  and  numerical  tech¬ 
niques  for  solving  large  body  coupling  and  scattering  problems.  The  second  involves  a ‘study  of  tech¬ 
niques  for  reducing  the  Radar  Cross-section  of  complex  targets  by  loading  them  with  either  a  resistive 
sheet,  or  with  a  lossy  dielectric  or  magnetic  material  that  may  be  arbitrarily  inhomogeneous.  The  third 
area  of  investigation  involves  the  study  of  radome  antennas  that  employ  frequency  selective  surfaces. 
Finally,  the  fourth  aspect  of  our  research  is  concerned  with  the  investigation  of  the  coupling,  cross-talk 
and  pulse  distortion  in  high-speed  digital  circuits. 

In  the  area  of  large  body  scattering  and  coupling  problems,  we  feel  that  we  have  made  significant 
progress  during  the  last  nine  months.  In  the  past,  we  have  been  involved  in  the  development  of  the 
spectral  domain  approaches  for  solving  electromagnetic  scattering  problems  above  the  resonance  region. 
During  the  course  of  this  research  it  was  determined  that  the  original  version  of  the  spectral  interac¬ 
tive  technique  (SIT),  while  it  is  very  versatile  and  numerically  efficient,  did  not  possess  the  feature  of 
guaranteed  convergence  for  arbitrary  structures.  This  was  not  totally  unexpected,  since  SIT  is  an 


asymptotic  technique,  similar  in  some  ways  to  the  geometrical  theory  of  diffraction,  or  the  GTD.  It  is 
well  known  that  the  asymptotic  methods  often  yield  a  divergent  series  when  the  higher  order  terms 
are  included.  However,  an  alternate  approach,  viz,  the  conjugate  gradient  (OG)  technique,  has  been 
recognized  for  its  salutary  feature  that  the  convergence  of  the  iteration  algorithm  employed  in  this  pro¬ 
cedure  is  guaranteed,  albeit  at  the  expense  of  slowing  down  the  iteration. 

During  the  last  nine  months  we  have  solved  a  large  class  of  representative  problems  by  combining 
the  moment  method,  as  well  as  the  spectral  domain  formulation,  with  the  conjugate  gradient  technique. 
We  found  that  while  for  many  cases  the  speed  of  convergence  the  conjugate  gradient  technique  was 
quite  adequate,  particularly  when  it  was  combined  with  the  spectra  domain  or  the  discrete  convolution 
techniques  that  allowed  the  use  of  the  FFT,  the  efficiency  had  considerable  room  for  further  improve¬ 
ment  for  some  large  body  coupling  and  scattering  problems.  In  light  of  this  fact,  we  proceeded  to 
develop  a  new  combinational  technique  that  attempts  to  combine  some  of  the  desirable  features  of  the 
SIT  and  the  CG  methods.  We  have  been  successful  in  doing  so,  and  we  are  now  able  to  guarantee  die 
convergence  of  the  spectral  iterative  technique  while  preserving  its  computational  efficiency. 

We  are  currently  in  the  process  of  testing  this  algorithm  by  applying  it  to  a  variety  of  different 
problems  involving  low  observable  targets  and  microwave  coupling  into  the  interior  of  vehicles,  eg, 
missiles. 

One  of  the  important  aspects  of  radome  antenna  design  is  to  lower  its  RCS  without  an  undue 
compromise  of  its  performance  as  an  antenna.  For  this  reason  we  have  studied  the  problem  of  FSS 
radomes  that  may  be  constructed  with  imperfectly  conducting,  e.g,  resistive  materials.  We  have 
developed  a  computer  program  for  studying  such  FSS  structures:  Other  studies  of  resistive  materials 
have  included  the  investigation  of  the  scattering  properties  of  bodies  that  are  made  up  of  resistive  cards, 
arranged  in  periodic  cells,  which  may  be  rectangular  or  hexagonal  in  shape.  The  attenuation  of 
electromagnetic  waves  propagating  through  such  lossy,  anisotropic  media  is  being  computed  together 
with  the  diffraction  properties  of  this  type  of  a  medium  as  seen  by  an  illuminating  radar. 

Finally,  we  have  investigated  the  crosstalk  and  pulse  distortion  in  multiconductor  transmission 
lines,  eg,  microstrip  lines  or  striplines,  when  these  lines  are  terminated  by  nonlinear  loads;  for  instance, 
logic  gates.  We  have  developed  a  new  time-domain  technique  for  solving  this  problem  which  cannot  be 
handled  using  the  conventional  frequency  domain  methods  that  employ  the  Fourier  transform  tech¬ 
nique,  since  the  nature  of  the  circuits  being  investigated  is  nonlinear.  We  are  currently  studying  mare 
complex  structures,  such  as  connectors  and  discontinuities  in  the  transmission  lines,  and  are  determining 
the  effects  of  these  discontinuities  on  the  degradation  of  high  speed  digital  pulses.  We  believe  that  the 
results  of  our  investigation  will  be  useful  in  the  problem  of  developing  CAD  tools  for  Packaging  of 
digital  circuits. 

SIGNIFICANT  RESEARCH  ACCOMPLISHMENTS: 

We  have  developed  a  hybrid  technique  based  upon  the  combination  of  the  Spectral  domain 
approach  and  the  minimization  of  an  error  functional  as  in  the  Conjugate  gradient  technique.  The  new 
algorithm  has  the  useful  feature  that  the  convergence  of  the  spectral  iteration  procedure  can  now  be 
guaranteed.  We  have  also  developed  a  technique  for  analyzing  inhomogeneously  coated  conducting 
bodies  and  bodies  that  are  constructed  from  resistive  cards.  A  new  time-supping  procedure  for  analyz¬ 
ing  digital  circuits  with  nonlinear  loads  has  also  been  developed. 

INTERACTIONS  AND/OR  TECHNOLOGY  TRANSFER: 

We  have  had  interactions  with  MICOM  at  Huntsville  who  have  been  interested  in  our  efforts  to 
build  a  220  GHz  imaging  radar.  We  have  transferred  the  design  of  such  a  radar  system  to  MICOM. 
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SCIENTIFIC  OBJECTIVE: 

Effective  utilization  of  the  millimeter-wave  (mm-wave)  spectrum,  for  communication,  radar  and 
electronic  counter  measure  and  for  several  other  applications,  depends  rather  critically  on  the  develop¬ 
ment  of  reliable,  low  cost,  planar  integrated  circuits.  A  major  objective  of  our  effort  is  the  development 
of  guiding  structures  and  associated  components  for  millimeter  and  submillimeter  wave  integrated  cir¬ 
cuits. 

Successful  development  of  these  waveguides  and  associated  components  requires  the  use  of  sophis¬ 
ticated  analytical  tools  capable  of  handling  complex  structures  comprising  dielectric  rods  and  substrates 
of  various  cross-section,  with  and  without  metallic  shields.  The  discontinuities  in  these  waveguides, 
that  are  inevitably  present  when  a  circuit  such  as  a  filter,  a  transition  or  an  oscillator  is  integrated  in  a 
monolithic  system,  pose  an  even  more  challenging  problem  to  the  designer  of  the  system.  The  availabil¬ 
ity  of  computer-aided  design  tools  offers  great  help  to  the  designer  and  a  major  objective  of  our  effort  is 
to  develop  such  tools. 

Since  few  reliable  and  efficient  analytical  techniques  are  currently  available  for  investigating 
planar  inhomogeneous  waveguides  and  discontinuities  that  can  often  support  higher  order  modes,  it  is 
essential  that  the  analytical  methods  we  develop  be  thoroughly  tested  experimentally.  For  this  reason, 
we  also  plan  to  pursue  an  experimental  program  that  would  support  and  complement  the  theoretical 
studies. 
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SUMMARY  OF  RESEARCH: 

During  the  past  nine  months  we  have  investigated  a  number  of  different  wave-guiding 
configurations  that  are  suitable  for  monolithic  millimeter  wave  circuits.  These  include  microstrip  lines, 
fin  lines  as  well  as  striplines  with  periodic  stubs.  The  design  of  monolithic  integrated  circuits  requires  as 
reliable  knowledge  of  the  equivalent  circuits  of  discontinuities  that  are  inevitably  present  in  these  cir¬ 
cuits.  With  this  in  mind,  we  have  been  investigating  a  number  of  different  techniques  for  analyzing 
and  characterizing  these  discontinuities  in  planar  waveguides.  We  have  also  been  pursuing  a  comple¬ 
mentary  experimental  program  to  verify  the  theoretical  results.  We  have  applied  the  finite  element 
technique  to  calculate  the  equivalent  circuit  of  discontinuities  in  planar  transmission  lines,  e.g,  a 
microstrip  line.  We  have  also  developed  the  mode-matching  solution  of  the  discontinuity  problem  in 
microstrips  and  fin  lines.  The  mode  matching  method  has  been  extensively  applied  in  the  past  to  solve 
the  discontinuity  problem  in  conventional,  hollow  metallic,  rectangular  or  circular  waveguides.  How¬ 
ever,  in  the  past,  a  lack  of  knowledge  of  the  higher  order  modes  in  planar  waveguides  prevented  the 
use  of  this  technique  for  the  solution  of  the  discontinuity  problem  in  a  microstrip  line  or  in  a  fin  line. 
For  this  reason,  we  have  developed  a  technique  for  accurate  computation  of  the  higher  order  modes  in 
such  waveguides.  The  availability  of  these  higher  order  modes  will  now  enable  us  to  solve  complex 
discontinuity  problems  using  the  mode  matching  procedure. 

In  addition  to  the  above  two  methods,  we  have  also  been  pursuing  the  development  of  the  spectral 
domain  approaches  for  solving  the  discontinuity  problems.  Iterative  techniques  are  currently  being 
developed  for  attacking  a  wide  class  of  discontinuity  problems  that  are  either  of  more  general  nature 
than  could  be  handled  with  the  mode  matching  technique,  or  are  too  time-consuming  or  unwieldy  to  be 
treatable  with  the  finite  element  method. 

Experimental  efforts  in  the  monolithic  integrated  circuit  area  have  included  the  modeling  of  the 
discontinuities  and  the  verification  of  the  theoretical  results;  Discontinuities  typically  present  in 
monolithic  designs,  e^»,  transitions  or  tapers  from  the  source  to  an  active  module,  or  from  the  source  to 
the  antenna,  are  being  investigated.  Other  types  of  discontinuities  being  studied  include  the  transition 
between  two  different  dielectric  substrates,  for  instance  duroid  (€,-2.2)  in  a  microstrip  line  and  the 
GaAs  (€,.-12.3)  substrate  in  the  active  module. 

In  addition  to  these  dielectric  antennas  compatible  with  the  planar  waveguides  have  been  experi¬ 
mentally  developed  for  operation  in  the  90  GHz  and  the  220  GHz  frequency  ranges.  These  antenna 
systems  have  been  used  to  measure  the  scattering  characteristics  of  a  number  of  canonical  geometries. 

In  the  area  of  quasi -optical  circuits  and  systems,  we  have  developed  a  technique  for  analyzing 
fiber  optical  guides  with  anisotropic  and  inhomogeneous  cores  that  have  special  single-mode  propagation 
characteristics.  We  have  also  investigated  mesh  filters  that  find  application  in  infrared  signature  work. 
These  filters  are  often  printed  on  substrates  whose  loss  characteristics  can  be  controlled  electrically, 
with  a  corresponding  alteration  of  their  reflection,  transmission  and  absorption  characteristics. 
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SCIENTIFIC  OBJECTIVE: 

Problems  of  control  systems  under  imperfect  information  arise  in  a  variety  of  applications.  Com¬ 
puter  control  of  aircraft,  electronic  guidance  and  navigation  of  missiles  and  submarines;  and  fire  control 
systems  are  typical  areas  of  application.  Uncertainties  arise  because  of  unknown  system  parameters, 
unknown  signal  environments;  and  hardware  tolerances.  Related  complexities  arise  in  situations 
involving  multiple  system  performance  criteria  and  multiple  decision  agents  having  access  to  different 
information  sets.  In  addition,  such  systems  may  consist  of  interconnections  of  systems,  Le,  be  modeled 
as  a  large  scale  system. 

The  principal  objective  of  this  research  unit  is  to  gain  basic  understanding  of  the  behavior  and 
control  of  complex  systems  containing  uncertainty.  Conceptual,  analytical,  and  computational  tools 
fundamental  to  the  synthesis  of  controllers  and  control  strategies  for  such  systems  will  be  developed. 
Stochastic  adaptive  systems,  robust  reduced  order  adaptive  systems,  sensitivity,  and  nonlinear  systems 
will  receive  specific  attention. 


SUMMARY  OF  RESEARCH: 

During  the  past  year,  we  have  continued  our  research  activities  on  the  fundamental  study  of  the 
role  of  modeling  and  information  uncertainties  in  complex  systems,  specifically  in  the  context  of  sto¬ 
chastic  adaptive  control,  stochastic  adaptive  Nash  strategies,  robust  reduced  order  adaptive  control,  sen¬ 
sitivity  and  robustness  of  multivariable  linear  systems,  and  nonlinear  systems. 

In  adaptive  control,  one  area  of  focus  has  been  the  development  of  a  new  stability  criterion  for 
systems  subjected  to  parameter  drifts.  Although  several  researchers  have  observed  parameter  drift 
instabilities  in  adaptive  schemes  for  identification  and  control,  an  analytical  explanation  of  this 
phenomenon  was  lacking.  We  have  recently  developed  a  theory  of  slow  adaptation  which  not  only 
explains  this  type  of  instability,  but  also  suggests  the  means  to  counteract  it.  Our  stability  criterion 
shows  an  explicit  dependence  on  input  signals  and  requires  that  most  of  signal  energy  be  in  the  fre¬ 
quency  range  where  the  real  part  of  the  transfer  function  is  positive.  The  criterion  has  attracted  the 


attention  of  many  researchers  not  only  in  this  country,  but  also  in  Australia  and  Sweden,  and  is  being 
extended  to  many  other  schemes  for  parameter  estimation  and  control.  In  another  study  in  the  adap¬ 
tive  field,  we  have  developed  an  optimization  oriented  approach  to  the  adaptive  control  of  systems 
described  by  Markov  Chains,  and  have  investigated  convergence  properties  of  the  resulting  algorithms. 
Furthermore,  we  have  extended  the  notion  of  adapation  to  multiple  decision  makers’  problems  under 
various  modes  of  decision  making,  including  Nash  and  hierarchical  solution  concepts.  Also,  the  Gen¬ 
eralized  Hessenberg  Representation  (GHR)  as  an  aggregation  oriented  model  of  large  systems  has  contin¬ 
ued  to  receive  development.  A  summary  of  recent  results  may  be  found  in  the  references  below.  Espe¬ 
cially  noteworthy  has  been  an  interpretation  via  the  GHR  of  multivariable  zeros  and  their  role  in  con¬ 
troller  design. 

During  this  period,  we  have  also  continued  work  on  the  theory  of  oscillators  and  on  ways  of 
using  the  decomposition  of  large  scale  systems  (both  controlled  and  uncontrolled)  to  deduce  various  pro¬ 
perties,  such  as  controllability  or  stability.  A  list  of  papers  describing  these  new  developments,  and 
published  either  in  the  second  half  of  1984  or  in  1985  (or  to  be  published)  is  provided  below. 
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SCIENTIFIC  OBJECTIVE: 

In  large  scale  systems,  control  tasks  may  be  decomposed  and  hierarchical  levels  may  be  imposed  on 
the  basis  of  analytically  derived  strategies.  This  project  is  devoted  to  a  study  of  such  hierarchical  stra¬ 
tegies,  and  to  associated  constraints  arising  from  computer  implementation  of  the  resulting  controllers. 
Questions  such  as  information  flow,  time  scales,  decentralization  and  loss  of  feedback  loops  will  be 
examined  in  coordinating  the  multiple  controllers  in  a  network. 


SUMMARY  OF  RESEARCH: 

During  the  past  year,  we  have  concentrated  our  efforts  primarily  on  (i)  time-scale  properties  of 
networks,  (ii)  integral  manifold  design  of  nonlinear  two-time  scale  systems,  (iii)  decentralized  control 
of  large  scale  systems  using  aggregation  techniques,  projective  controls  or,  dominance,  (iv)  computer- 
aided-design  of  control  systems,  (v)  minima*  design  philosophy  and  its  applications  to  stochastic  control 
systems  with  unknown  parameters,  and  to  communication  and  jamming  problems  with  partially 
unknown  statistics,  and  (vi)  developing  decentralized  algorithms  for  multiple  decision  maker  problems. 

For  graphs  representing  massive  nodes  connected  with  0-1  branches,  we  have  developed  a 
methodology,  showing  how  to  represent  them  in  a  simpler  aggregate  form.  This  methodology  exploits  a 
newly  established  fact  about  the  relationship  between  the  connection  density  and  graph  eigenvalues. 
Using  the  eigenvalues  as  time-scale  indicators,  the  aggregate  graph  interpretation  is  that  it  represents 
slow  system-wide  phenomena,  while  the  faster  phenomena  are  localized  in  the  areas  of  denser 
connections.  Several  tests  on  graphs  of  large  scale  systems  confirm  this  modeling  procedure. 

In  the  field  of  nonlinear  systems,  by  approaching  the  design  of  two-time  scale  systems  in  two 
stages,  first,  to  reach  the  manifold,  and,  second,  to  satisfy  a  performance  requirement  in  the  manifold. 
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we  have  been  able  to  solve  several  nonlinear  problems  of  practical  interest  Among  them  are  the  feed¬ 
back  linearization  problem  for  a  flexible  robotic  manipulator,  and  a  design  method  fo^.  systems  with 
imperfect  actuators.  For  a  manipulator  with  flexible  joints,  our  design  achieves  a  performance  close  to 
that  of  a  perfectly  rigid  manipulator.  These  designs  have  validated  the  general  approach  suggested  by 
one  of  us  (PVK)  during  a  plenary  talk  at  the  1984  IF  AC  World  Congress,  that  the  concept  of  an 
integral  manifold  can  be  used  to  unify  most  nonlinear  studies. 

On  the  topic  of  minimax  design  philosophy,  we  have  developed  a  number  of  theoretical  results 
verifying  existence  of  saddle  points  for  stochastic  optimization  problems  with  nonclassical  decentralized 
information,  with  one  source  of  motivation  being  jamming  problems  arising  in  communications.  We 
have  also  applied  this  philosophy  to  the  design  of  controllers  for  stochastic  systems  with  a  hybrid  type 
of  uncertainty,  and  have  conducted  a  numerical  study  on  the  comparison  between  performance  meas¬ 
ures  under  the  minimax  and  the  traditional  stochastic  design  techniques. 

Furthermore,  development  of  the  computer-aided  design  package  L-A-S  continues.  L-A-S  pro¬ 
vides  both  design  support  for  simulation  experiments  with  large  scale  systems,  and  a  forum  for 
algorithm-oriented  control  law  synthesis  as  well. 

A  list  of  papers  describing  these  or  earlier  research  results,  and  published  either  in  1984  or  in  1985 
(or  to  be  published)  is  provided  below. 
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SCIENTIFIC  OBJECTIVE: 

The  objectives  of  the  proposed  research  are  to  develop  new  analytical  methods  for  performance 
evaluation  of  multiple-terminal  digital  communication  systems;  devise  new  techniques  for  modulation, 
coding,  and  diversity  which  will  provide  increased  jamming  protection  and  multiple-access  capability 
in  spread-spectrum  systems;  evaluate  alternative  receiver  structures  and  diversity  combining  methods 
for  spread-spectrum  communication  in  a  fading  environment;  and  investigate  new  approaches  to 
transmission  scheduling  and  routing  in  mobile  radio  networks. 


SUMMARY  OF  RESEARCH: 

/.  Diversity  Combining  for  Frequency-Hop  Communications 

Diversity  transmission  is  often  employed  in  frequency-hop  systems  to  provide  reliable  communi¬ 
cation  in  the  presence  of  fading  or  partial-band  interference.  Some  of  the  methods  for  combining  the 
diversity  receptions  require  side  information  regarding  the  presence  or  absence  of  interference  on  the 
diversity  receptions.  Because  of  implementation  considerations,  diversity  combining  techniques  that  do 
not  require  side  information  from  external  sources  are  particularly  attractive.  One  diversity  combining 
scheme  that  does  not  require  side  information  is  called  dipped  linear  combining  [l  J.  In  this  system,  the 
output  of  the  envelope  detector  of  each  diversity  reception  is  clipped  before  it  is  added  with  the  other 
diversity  receptions.  It  is  shown  in  [l]  that  clipped  linear  combining  is  effective  against  partial-band 
interference.  There  is  a  practical  disadvantage  to  this  diversity  combining  scheme,  however,  since  the 
clipping  level  depends  on  the  average  received  signal  voltage  in  the  absence  of  the  interference.  This 
signal  voltage  may  be  difficult  to  measure  in  practice. 

It  is  desirable  to  employ  diversity  combining  techniques  that  do  not  depend  on  the  received  signal 
power.  One  such  diversity  combining  technique  is  ViterWs  ratio  threshold  technique  [2}-[4j.  Each 
diversity  reception  is  determined  to  be  accepted  or  rejected  based  on  the  relative  values  of  the  envelope 
detector  outputs.  More  specifically,  for  each  diversity  reception,  the  ratio  of  the  largest  envelope  detec¬ 
tor  output  to  the  second  largest  envelope  detector  output  is  compared  to  a  threshold  which  is  a  fixed 
number  greater  than  one.  If  the  ratio  is  less  than  the  threshold,  the  diversity  reception  is  rejected,  and 
if  it  is  greater  than  the  threshold,  it  is  accepted.  If  at  least  one  reception  is  accepted,  then  only  the 
accepted  diversity  receptions  are  combined.  If  all  diversity  receptions  are  rejected,  then  all  of  them  are 
combined  and  a  hard  decision  is  made  (the  resulting  symbol  may  be  flagged  for  erasure  if  the  decoder 
has  a  provision  for  correcting  erasures  [5]X 

We  have  computed  the  average  probability  of  error  for  Viterbi’s  ratio  threshold  technique  used  in 
conjunction  with  diversity  for  frequency-hop  communications  in  the  presence  of  partial-band  Gaussian 
interference.  The  system  with  binary  orthogonal  signaling  and  noncoherent  demodulation  is  analyzed. 
A  comparison  is  made  between  the  system  employing  Viterbi’s  ratio  threshold  technique  with  linear 


diversity  combining  and  the  system  that  utilizes  clipped  linear  combining  [l]  as  the  diversity  scheme. 
We  allow  a  nonzero  background  or  quiescent  noise  level  to  account  for  thermal  noise  in  the  receiver 
and  other  wideband  noise  sources. 

2  Fast  Frequency  Hopping  for  Mobile  Radio  Applications 

We  have  continued  out  investigation  of  the  performance  of  a  class  of  fast  frequency  hopped 
spread-spectrum  multiple-access  communication  schemes  proposed  for  mobile  radio  communications. 
We  have  studied  the  performance  of  the  sub-optimum  receiver  (derived  in  earlier  work)  when  it  uses 
hard  decisions  instead  of  soft  decisions  on  each  hop.  Both  the  additive  white  Gaussian  noise  channel 
and  the  nonselective  Rayleigh  fading  channel  have  been  considered.  We  have  also  analyzed  the  perfor¬ 
mance  of  this  receiver  when  interfering  signals  from  other  users  are  present  Analytical  results  have 
been  obtained  for  both  the  nonfaded  and  the  nonselective  Rayleigh  faded  channel.  These  results  indi¬ 
cate  that  the  performance  of  this  receiver  is  somewhat  different  than  the  performance  of  previously 
considered  receivers  for  this  system  [6]. 

3.  Spread-Spectrum  Multiple- Access  Communication  for  the  Low  Frequency  Band 

We  have  continued  our  investigation  of  DS/SSMA  communication  over  the  low  frequency  (LF) 
band.  In  previous  work,  we  had  derived  a  locally  optimum  receiver  for  such  a  system  operating  in  a 
mixture  of  Gaussian  and  impulsive  noise.  The  receiver  consists  of  a  bandpass  correlator  followed  by  a 
sampler,  a  zero-memory  nonlinearity,  and  M  discrete-time  matched  filters.  More  recently,  we  have  stu¬ 
died  the  performance  of  such  a  receiver  when  the  optimum  nonlinearity  is  replaced  by  mare  easily 
realizable  nonlinearities  such  as  limiting  amplifiers  hard  limiters  and  hole-punchersL  Our  results  indi¬ 
cate  that  performance  (as  measured  by  the  signal-to-noise  ratio)  is  degraded  only  slightly  when  such 
nonlinearities  are  used,  and  that  the  performance  is  not  very  critically  dependent  on  the  parameters 
defining  these  nonlinearities.  Error  probability  as  a  performance  measure  is  also  being  considered. 

4.  Transmission  Scheduling  and  Retransmission  Control  in  Spread-Spectrum  Radio  Networks 

We  have  made  more  progress  on  the  problem  of  scheduling  transmissions  in  a  spread-spectrum 
mobile  packet  radio  network  to  avoid  primary  collisions.  We  have  recently  discovered  a  polynomial 
time  algorithm  for  computing  a  schedule  of  minimal  length.  This  minimal  length  is  also  the  mul- 
tichromatic  index  of  a  graph.  The  algorithm  is  based  on  work  of  Grotschel,  Lovasz,  and  Schrijver  and 
work  of  Edmonds.  The  algorithm  proves  false  the  previously  held  belief  that  the  scheduling  problem 
is  NP-complete. 

In  another  line  of  JSEP  funded  research,  we  have  shown  how  the  theory  of  stochastic  approxima¬ 
tion  can  be  used  to  fruitfully  study  adaptive  control  of  retransmission  probabilities;  A  quite  general 
family  of  retransmission  control  strategies  was  introduced  such  that,  on  the  basis  of  very  limited  feed¬ 
back,  the  transmission  probabilities  of  various  stations  can  be  controlled. 

Issues  of  fairness  and  priorities  are  addressed  through  the  introduction  of  a  cost  function  which 
the  stations  endeavor  to  minimize  in  a  distributed  fashion.  The  estimators  appear  to  be  quite  effective 
in  simulations.  The  use  of  stochastic  approximation  provides  a  fairly  systematic  approach  which  can  be 
applied  as  well  to  the  control  of  transmission  power,  code  rates,  etc.  These  possibilities  deserve  further 
study. 

5.  Clustering  Algorithms  for  Packet  Radio  Networks 

We  have  begun  a  study  of  algorithms  for  organizing  a  packet  radio  network  into  overlapping 
clusters.  Two  strategies  were  proposed  and  compared.  The  strategies  incorporate  heuristics  for  generat¬ 
ing  priorities  which  determine  which  nodes  will  become  cluster  heads.  Simulations  were  used  to  gen¬ 
erate  trade-off  curves  for  the  strategies.  A  small  mean  distance  to  the  nearest  cluster  head  can  be  traded 
off  to  obtain  a  low  density  of  cluster  heads.  As  we  found,  some  clustering  strategies  tend  to  have  more 
favorable  trade-off  curves. 
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SCIENTIFIC  OBJECTIVES: 

The  primary  goal  of  this  research  is  the  development  of  practical  and  efficient  implementations 
for  signal  detection  and  estimation  systems  with  potential  applications  to  communication  receivers, 
sonar  and  radar  tracking  and  detection,  and  other  statistical  signal  processing  functions.  The  main 
emphasis  is  on  the  development  of  techniques  for  the  design  and  analysis  of  digital  systems  for  use  in 
such  applications,  particular  attention  will  be  paid  to  situations  in  which  signal  and  noise  models  do 
not  obey  standard  assumptions  or  in  which  accurate  modeling  is  not  possible.  Such  situations  invite  the 
utilization  of  nonlinear,  robust,  adaptive,  and/or  nonparametric  systems  which,  while  exhibiting  supe¬ 
rior  performance,  often  require  complex  structures  for  exact  implementation.  Advances  in  integrated 
circuitry  have  made  practical  the  approximation  of  such  structures  with  digital  systems,  and  thus  this 
study  is  motivated  in  part  by  the  desire  to  utilize  fully  the  flexibility  and  power  of  the  new  technol¬ 
ogy. 


SUMMARY  OF  RESEARCH: 

Our  JSEP  sponsored  research  during  this  reporting  period  has  been  focused  primarily  on  two  prob¬ 
lems:  (l)  performance  evaluation  of  direct-sequence  spread-spectrum  multiple-access  (DS/SSMA)  com¬ 
munications  in  non-Gaussian  noise  environments;  and  (2)  analysis  of  data  quantization  schemes  in  sig¬ 
nal  detection  and  estimation.  Also,  during  this  period  we  have  prepared  an  extensive  survey  for  the 
Proceedings  of  the  IEEE  describing  our  work  on  the  problem  of  robust  filtering  that  has  been  con¬ 
ducted  under  JSEP  and  other  sponsorship  over  the  past  six  years. 

Non-Gaussian,  impulsive  noise  is  a  key  noise  source  in  many  multiple-access  communication  sys¬ 
tems  (such  as  military  radio  networks).  Thus,  when  studying  such  systems,  the  evaluation  of  their  per¬ 
formance  in  non-Gaussian  noise  is  of  interest.  We  have  considered  two  techniques  for  evaluating  the 
average  error  probabilities  of  conventional  linear  DS/SSMA  receivers  in  impulsive  noise,  one  of  which 
uses  a  characteristic  function  method  and  is  useful  for  small  spreading  sequence  lengths  and  the  other 
of  which  uses  a  Taylor-series  expansion  method  for  longer  spreading  sequences.  We  carry  these  results 
further  by  considering  nonlinear  correlation  receivers  to  improve  the  receiver  performance  against 
impulsive  noise.  These  results  indicate  that  substantial  performance  gain  can  be  obtained  by  using  non¬ 
linear  reception,  particularly  for  long  spreading  sequences  and  heavily  impulsive  noise. 

We  also  have  considered  several  problems  relating  to  the  quantization  of  data  for  various  pur¬ 
poses.  In  this  work  we  have  concentrated  on  the  problem  of  fine  quantization  and  on  the  development 
of  performance  expressions  and  design  methods  for  this  case.  Specific  problems  that  we  have  treated  are 
the  analysis  and  design  of  quantizers  for  signal-detection  data,  the  analysis  of  the  effects  of  data  quanti¬ 
zation  on  MMSE  filtering,  and  the  design  of  quantizers  that  have  good  distortion  performance  for  data 
with  uncertain  statistics. 
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A  further  study  has  treated  the  analysis  of  the  efficiencies  of  a  class  of  multistage  signal-detection 
algorithms.  These  algorithms  are  seen  in  our  analysis  to  have  efficiencies  near  that  of  the  optimum 
sequential-probability-ratio  detector  while  retaining  the  design  and  implementational  simplicity  of 
conventional  fixed-sample-size  detection  algorithms.  These  new  algorithms  are  seen  to  be  particularly 
efficient  for  applications,  such  as  search  radar,  in  which  the  presence  of  a  signal  occurs  rarely  relative  to 
the  absence  of  one. 
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SCIENTIFIC  OBJECTIVE: 

Since  the  introduction  of  the  integrated  circuit  chip,  chip  complexity  has  doubled  every  two 
years.  However,  concurrent  with  this  growth  in  complexity,  the  design  cycle  has  lengthened  and  the 
design  expense  has  been  increasing  by  a  factor  of  three  or  more  about  every  five  years:  We  are  now  at 
the  point  where  the  design  of  a  custom  VLSI  circuit  will  require  over  one  hundred  man-years  with 
present  design  methodologies.  Thus,  the  initiation  of  a  new  VLSI  circuit  is  a  very  major  investment 
risk.  This  trend  in  design  complexity  is  a  major  concern  of  leaders  in  the  electronics  industry.  New 
design  methodologies  as  well  as  better  computer-aids  for  the  synthesis,  layout,  simulation,  verification, 
and  testing  of  VLSI  systems  are  needed.  Unfortunately,  the  computational  complexity  of  the  problems 
in  this  area  are  typically  NP-complete.  From  the  point  of  view  of  the  development  of  design  automa¬ 
tion  tools,  the  apparent  intractability  of  NP-complete  problems  leads  one  to  the  use  of  suboptimal 
heuristic  algorithms.  Because  of  this  complexity  a  hierarchical  approach  will  have  to  be  taken  in  the 
solution  of  the  problems  in  this  field. 

The  objective  of  our  research  in  this  area  is  to  study  hierarchical  approaches  to  the  simulation  and 
verification  of  VLSI  circuit  designs.  In  particular,  models  and  design  verification  aids  for  the  various 
hierarchical  levels  will  be  investigated.  The  trend  in  VLSI  design  is  to  use  a  top-down  structured 
design  approach  in  which  systems  are  implemented  by  means  of  a  number  of  identical  cells  which  com¬ 
municate  to  each  other  in  a  predefined  way.  It  is  important  in  the  bottom-up  design  verification  phase 
to  be  able  to  model  and  accurately  simulate  the  timing  of  these  communication  links,  particularly  in 
high  speed  signal  processing  systems. 
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SUMMARY  OF  RESEARCH: 

This  past  year  a  new  DC  circuit  model  for  modulation  doped  GaAs  FETs  was  developed  [1,10]. 
Other  models  are  unnecessarily  complex  and  have  discontinuous  derivatives  at  the  boundaries  between 
the  different  regions  of  operation.  Thus,  in  circuit  simulation  the  models  are  very  costly.  Furthermore, 
past  experience  with  MOSFET  models  has  shown  that  models  with  discontinuities  in  their  derivatives 
have  serious  convergences  problems.  Unfortunately,  only  recently  have  researchers  discovered  that 
these  discontinuities  were  the  source  of  this  very  serious  problem  which  has  plagued  simulation  experts 
for  a  decade.  Our  new  model  overcomes  these  difficulties,  and  it  is  simple  and  accurate.  The  model  can 
be  used  for  both  MODFETs  (HEMTs)  or  MOSFETs,  and  even  though  the  model  parameters  have 
different  interpretations  for  each  of  these  transistor  types,  physical  arguments  can  be  made  as  to  why 
the  model  can  be  used  for  both  devices  [l]. 

The  model  has  been  implemented  in  a  parameter  characterization  system  which  consists  of  a  HP 
4145  Semiconductor  Parameter  Analyzer  and  a  HP  9836  computer.  An  optimization  routine  is  used  to 
fit  the  model  parameters  to  the  measured  device  characteristics.  Excellent  agreement  has  been  obtained 
between  the  measured  characteristics  and  the  model  on  several  MODFET  samples  provided  by  Professor 
Morkoc. 

Work  is  continuing  on  the  hierarchical  simulation  of  VLSI  circuits  for  timing  analysis.  Several 
years  ago  under  JSEP  sponsorship  we  used  node  tearing  methods  to  partition  large  circuits  hierarchi¬ 
cally.  This  allowed  us  to  (a)  use  block  LU  factorization  techniques  in  the  solution  of  the  circuit  equa¬ 
tions,  (b)  declare  circuits  latent  at  higher  levels  in  the  hierarchy,  and  (c)  reduce  memory  requirements. 
The  block  LU  factorization  of  each  subcircuit  equations  can  proceed  in  parallel  on  a  multiprocessor  sys¬ 
tem.  Recently  a  study  was  completed  on  the  use  of  special  purpose  computer  architectures  for  the  LU 
factorization  of  this  partitioned  system  of  equations  [ill  The  performance  of  both  systolic  and  wave- 
front  array  processors  was  investigated.  However,  the  LU  factorization  of  the  interconnect  equations 
remains  a  bottleneck  in  the  effort  to  obtain  significant  speed  improvement.  Thus,  our  attention  turned 
to  decoupling  methods.  Unfortunately  these  methods  add  another  iteration  loop  to  the  solution  process, 
and  one  must  be  concerned  about  convergence  within  this  loop.  We  studied  the  convergence  properties 
of  two  decoupling  methods,  an  incremental  Gauss  Seidel  relaxation  algorithm  and  the  waveform  relax¬ 
ation  algorithm.  In  the  algorithms  used  in  standard  circuit  simulators  the  convergence  depends  on  the 
local  truncation  error  and  thus  the  timestep  can  be  controlled  by  the  rate  of  change  of  the  waveform. 
However,  when  relaxation  techniques  are  employed,  we  found  that  the  convergence  is  now  a  function 
of  the  coupling  among  the  partitioned  circuits.  If  the  coupling  is  sufficiently  weak  the  time  step  can 
still  be  controlled  by  the  local  truncation  error.  However,  it  is  not  always  easy  to  determine  the  degree 
of  coupling  and  so  the  time  step  must  not  only  be  controlled  by  the  local  truncation  error  but  also  by 
the  number  of  iterations  in  order  to  achieve  good  convergence  [2,8,9],  In  cases  where  strong  coupling 
exists  the  method  can  be  extremely  costly.  Thus,  the  performance  of  these  relaxation  algorithms  is  crit¬ 
ically  dependent  on  the  partitioning  of  the  circuit. 

Earlier  it  was  pointed  out  that  in  the  solution  of  nonlinear  circuit  equations,  the  convergence  of 
the  algorithm  depends  on  the  functional  description  of  the  model.  We  have  found  that  piecewise  linear 
solution  techniques  result  in  very  robust  algorithms  [4,5].  These  techniques  require  that  the  nonlinear 
device  characteristics  be  described  by  a  series  of  piecewise  linear  segments.  This  model  is  more  flexible 
and  more  readily  adaptable  to  changing  technology,  however,  the  description  is  not  in  terms  of  physical 
parameters.  It  may  be  necessary  to  use  a  physical  model  description  externally,  but  internally  a  piece- 
wise  linear  description  could  be  generated  in  order  to  achieve  reliable  convergence  to  the  solution. 

Finally,  several  publications  have  appeared  on  work  in  the  design  verification  area.  This  work 
was  initially  funded  by  JSEP.  The  goal  was  to  develop  algorithms  to  automatically  partition  a  transis¬ 
tor  netlist  description  of  a  VLSI  circuit  into  functional  blocks,  to  then  generate  functional  expressions 
to  describe  the  logic  behavior  of  each  block  [3l  and  to  generate  timing  models  [12].  These  algorithms 
have  been  implemented  into  two  programs  EXPRESS  and  MOSTIM  under  industrial  support. 

In  order  to  generate  accurate  timing  information  it  is  essential  to  have  good  circuit  models  for  the 
interconnect.  Thus,  under  JSEP  support  work  was  initiated  on  the  generation  of  RC  models  for  the 
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interconnect.  Recently  under  industrial  support  some  of  these  models  were  compared  to  measured 
results  from  several  test  wafers  [l\  Taking  into  account  tolerance  effects,  the  measured  results  for  both 
se  If  -capacitance  and  coupling  capacitance  of  the  interconnect  compared  favorably  with  the  numbers 
predicted  by  the  model. 


INTERACTIONS  AND/OR  TECHNOLOGY  TRANSFER: 

We  have  5  to  10  request  a  year  from  industry  and  universities  for  the  circuit  simulation  program 
SLATE  SLATE  has  the  same  input  format  as  SPICE  but  is  a  hierarchical  simulator  which  typically 
runs  2  to  3  times  faster  than  SPICE  and  requires  half  the  memory.  It  also  has  much  better  convergence 
than  SPICE2  because  of  an  improved  MOSFfcT  model 

As  a  result  of  our  research  on  interconnect  models,  Intel  Corporation  is  finding  our  results  useful 
to  estimate  interconnect  capacitance  in  VLSI  circuits. 


PUBLICATIONS: 


Journal  Articles: 

1.  KJl  Cioffi  and  T.N.  Trick,  "A  DC  circuit  model  far  modulation  doped  FETV  submitted  to 
IEEE  Electron  Device  Letters.  (JSEP) 


Conference  Papers : 

2.  WJC  Chia,  TJM.  Trick,  and  LN.  Haji  "Implementation  of  a  relaxation  technique  in  a  general  pur¬ 
pose  circuit  simulator,"  Proc.  of  Int.  Conference  on  Circuits  and  Systems,  June  10-12, 1985,  to 
appear.  (JSEP) 

1  LN.  Hajj  and  D.  Saab,  "On  the  function  logic  representation  of  digital  transistor  circuits,"  IEEE 
Int.  Symp.  on  Circuits  and  System,  Kyoto,  Japan,  June  1985.  (SRC/JSEP) 

4.  LG.  Mao  and  LN.  Hajj,  "A  piecewise- linear  approach  to  dc  analysis  of  large-scale  integrated  cir¬ 
cuits,"  The  Fifth  Int.  Symp.  on  Network  Theory,  Sarajevo,  Yugoslavia,  pp-  472-476,  September 
1984.  (JSEP) 

5.  LG.  Mao  and  LN.  Haji  "A  Gauss-Seidel  piecewise- linear  method  for  the  analysis  of  large-scale 
integrated  circuits,"  International  Conf.  on  Circuits  and  Systems,  Beijing,  China,  June  1985. 
(JSEP) 

6.  DuC.  Yeh  and  TJM.  Trick,  "Suboptimal  and  optimal  algorithms  for  PLA  folding,"  Proc.  IEEE 
Symp.  on  Circuits  and  systems,  Kyoto,  Japan,  June  5-7, 1985,  to  appear.  (JSEP) 

7.  CP.  Yuan  and  TJM.  Trick,  "Calculation  of  capacitance  in  VLSI  circuits,"  Proc.  IEEE  Int.  Conf. 
on  Computer-Aided  Design,  pp.  263-265,  Nov.  12-15, 1984.  (JSEP/IBM) 

8.  TJM.  Trick,  (Invited  talk)  "Experiences  using  a  relaxation  based  simulator  for  timing  analysis,” 
Spring  ACM  Design  Automation  Workshop,  January  21-22, 1985.  (JSEP/SRC) 


Theses: 


9. 


WJC  Chia,  "A  study  of  relaxation  techniques  for  the  transient  analysis  of  digital  circuits,"  PhD. 
Thesis,  August  1984.  (JSEP/SRC) 
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10.  KJL  Cioffi,  "Modeling  modulation  doped  field  effect  transistors  for  use  with  the  TECAP  systems," 
MS.  Thesis,  May  1984.  (JSEP) 

11.  V.B.  Rao,  "Switch-level  timing  simulation  of  MOS  VLSI  circuits,"  Phi).  Thesis,  January  1985. 
(JSEP/IBM) 

12.  K.-M.  Ivy  Lui,  "Special  purpose  computer  architecture  for  LU  factorization  of  partitioned  sys¬ 
tems,"  MS.  Thesis,  August  1984.  (JSEP) 


PUBLICATIONS  (OTHER  SUPPORT); 

Journal  Articles 

1.  D.E.  Hocevar,  MIL  Lightner  and  TAT.  Trick,  "An  extrapolated  yield  approximation  technique  for 
yield  maximization,"  IEEE  Trans,  on  Computer-Aided  of  Integrated  Circuits  and  Systems, 
CAD-3,  pp.  279-287,  Oct.  1984.  (NSF) 

2.  DJL  Hocevar,  P.  Yang,  TJM.  Trick,  and  BJP.  Epler,  "Transient  sensitivity  computation  for  M05FET 
circuits,"  IEEE  Trans,  on  Computer-Aided  Design  of  Integrated  Circuits  and  Systems,  CAD-4, 
October  1985,  to  appear. 

Conference  papers 

1.  LN.  Haj>  D.  Saab  and  B.  Rosario,  "Logic  and  timing  simulation  of  bipolar  ECL  circuits,"  1984  Int. 
Conf.  on  Computer-Aided  Design,  Santa  Clara,  CA,  pp.  194-196,  November  1984.  (SRC,  TRW, 
TEKTRONIX) 

2.  VS  Rao  and  TJSI.  Trick,  "A  new  approach  to  processing  strongly  connected  circuit  blocks  in  a 
waveform  relaxation  switch-level  timing  simulator,"  Proc.  IEEE  Int.  Conf.  on  Computer 
Design:  VLSI  in  Computer,  pp.  502-506,  Oct.  8-11, 1984.  (IBM) 

3.  V3.  Rao  and  T.N.  Trick,  "Switch-level  simulation  of  MOS  VLSI  circuits,"  Proc.  IEEE  Symp.  on 
Circuits  and  Systems,  June  5-7, 1985,  to  appear.  (IBM). 

4.  D.  Saab  and  LN.  Hajj,  "Parallel  and  concurrent  fault  simulation  of  MOS  circuits,"  IEEE  Int. 
Conference  on  Computer  Design  ( ICDD84),  Port  Chester,  NY,  pp.  752-756,  October  1984.  (SRC) 
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Fellow,  IEEE 
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WORK  UNIT  NUMBER  20 


TITLE:  Computer  Architecture 


SENIOR  PRINCIPAL  INVESTIGATORS: 

Edward  &  Davidson,  Research  Professor 
Janak  H.  Patel,  Research  Associate  Professor 
Prithviraj  Banerjee,  Research  Assistant  Professor 


SCIENTIFIC  PERSONNEL  AND  TITLES: 

Hoichi  Cheong,  Research  Assistant 
Daniel  J.  Colg  lazier.  Research  Assistant 
Richard  J.  Eickemeyer,  Research  Assistant 
Peter  Y.-T.  Hsu,  Research  Assistant 
Subhasis  Laha,  Research  Assistant 
Mohammad  Malkawi,  Research  Assistant 
Geoffrey  D.  McNiven,  Research  Assistant 
Gurindar  &  So  hi.  Research  Assistant 
Shu-Hui  Sung,  Research  Assistant 


SCIENTIFIC  OBJECTIVES: 

This  unit  seeks  to  develop,  model,  and  analyze  efficient,  high-performance  computer  architectures. 
We  have  identified  several  aspects  of  computer  architecture  for  particular  attention  due  both  to  their 
emerging  importance  from  a  technology-driven  point  of  view  and  to  the  lack  of  known  structures  or 
analysis  techniques  for  meeting  our  objective. 

1.1  Parallel- pi  peline  architectures 

This  section  is  devoted  to  exploring  the  limits  of  highly  concurrent  VLSI  processors  and  multipro¬ 
cessors.  VLSI  technology  affords  large  on-chip  complexity  but  has  relatively  limited  capability  for 
interchip  communication.  Accordingly,  there  is  a  heavy  emphasis  on  using  on-chip  storage  effectively 
to  reduce  off-chip  communication.  Furthermore,  the  long  access  times,  albeit  with  high  bandwidth, 
associated  with  off-chip  requests  creates  the  need  for  issuing  such  requests  earlier  to  achieve  effective 
forwarding  of  information  to  a  processor  so  as  to  minimize  the  wait  time  of  a  highly  pipelined  proces¬ 
sor.  Finally,  dependencies  and  linkages  between  tasks  allocated  to  different  processors  must  be  studied 
to  achieve  minimal  dependency-related  performance  degradation. 

13.  Memory  organizations 

In  recent  years,  the  memory  has  become  one  of  the  most  expensive  parts  of  a  typical  computing 
system.  In  addition,  the  overall  performance  of  a  system  is  also  primarily  governed  by  its  memory 
organization.  Therefore,  the  investigation  of  memory  organizations  is  essential  to  the  development  of 
high-performance  computing  systems.  In  particular,  multiprocessor  systems  are  becoming  increasingly 
important  due  to  the  cost-effectiveness  considerations  of  VLSI  technology.  To  exploit  this  potential. 
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memory  system  architectures  must  be  appropriately  organized  for  such  systems.  The  objective  of  this 
research  is  to  develop,  model,  and  analyze  high-performance  memory  organizations  for  multiprocessing 
and  multiprogramming  systems,  including  cache  and  virtual  memories,  for  symbolic  and  numeric  pro¬ 
cessing. 

SUMMARY  OF  RESEARCH: 

1.1  Parallel-pipeline  architectures 

We  have  developed  a  program-referencing  analysis  facility  that  determines  the  best  use  of  an 
expanded  capability  for  on-chip  storage  for  the  purpose  of  minimising  off-chip  memory  references. 
Optimal  tradeoffs  between  storing  instructions,  memory  accessing  parameters,  program  control  parame¬ 
ters,  and  operands  can  be  found  and  the  effectiveness  of  various  amounts  of  on-chip  storage  can  be 
determined. 

We  have  developed  a  formal  model  that  characterizes  the  information-accessing  needs  of  program 
modules  at  various  levels  of  granularity.  State-saving  requirements  for  rapid  exception  and  error 
recovery,  as  well  as  insights  for  scheduling  and  allocation  of  modules  to  processors  and  effective 
information-forwarding  techniques  will  be  derived  from  the  use  of  this  model. 

A  new  architecture  for  LISP  machines  is  being  developed  that  employs  small  associative  tables  for 
storing  S-expressions.  This  representation  is  compact  and  breaks  through  the  traditional  serial-execution 
barrier  for  LISP  posed  by  the  traditional  representation  using  indirect  pointer  chains.  Parallel  garbage 
collection  and  information  forwarding  capabilities  will  allow  highly  concurrent  pipelined  LISP 
machines  to  be  developed. 

1.2  Memory  organizations 

We  have  studied  several  numerical  programs  in  the  context  of  our  research  on  virtual  memory 
management.  Our  study  revealed  that  program  localities  can  be  quantitatively  reconized  through  some 
analysis  of  source  code.  Localities  so  identified  can  then  be  made  into  memory  management  directives 
to  be  used  at  execution  time.  This  is  the  basic  idea  behind  the  compiler-directed  (CD)  memory  manage¬ 
ment  policy.  Our  initial  experiments  in  a  uniprocessor,  uniprogramming  environment  show  that  the 
CD  policy  out-performs  the  well-known  least-recently  used  (LRU)  replacement  policy  and  the  -work¬ 
ing  set  policy  in  terms  of  number  of  page  faults  and  dynamic  memory  usage. 

For  list  processing,  we  have  developed  a  parallel  garbage  collection  scheme  and  the  underlying 
architecture  which  avoids  the  synchronization  overhead  found  in  most  other  parallel  schemes.  More¬ 
over,  in  our  method,  the  page  faults  caused  by  the  garbage  collector  do  not  adversely  affect  the  list  pro¬ 
cessor  performance,  unlike  in  other  schemes. 


PUBLICATIONS: 

Journal  Articles: 

1.  T.C.-K.  Chou  and  J.A.  Abraham,  "Distributed  control  of  computer  systems,"  IEEE  Trans.  Corn- 

put.  { to  appear).  (JSEP/VHSIC) 

Conference  Papers: 

1.  P.Y.-T.  Hsu,  J.T.  Rahmeh,  ES.  Davidson,  and  LA.  Abraham,  "TIDBITS:  Speedup  via  time-delay 

bit-slicing  in  ALU  design  for  VLSI  technology,  Proc.  12th  Int.  Symp .  Com  put.  Arch,  June  17-19, 
1985  (to  appearX  (SRC/JSEP) 


2.  A.  Ram  and  JJL  Patel,  "Parallel  garbage  collection  without  synchronization  overhead,"  Proc.  12th 
Int.  Symp.  Com  put.  Arc.,  June  17-19, 1985  (to  appear).  (JSEP) 

3.  G.S.  Sohi,  ES  Davidson,  and  JJL  Patel,  "An  efficient  LISP-execution  architecture  with  a  new 
representation  for  list  structures,"  Proc.  12th  Ira.  symp.  Comput.  Arch*  June  17-19,  1985  (to 
appear).  (JSEP) 

4.  WA  Sufah,  R.  Lee  and  M.  Malkawi,  "Identifying  two  program  categories  for  memory  manage¬ 
ment  purposes,"  Proc.  8th  Int.  Conf.  Comput.  Software  and  Applications,  Nov.  7-9,  1984. 
(NSF/DOE/JSEP) 


Technical  Reports  and  Documents: 


1.  D.  Colglazier,  "A  performance  analysis  of  multiprocessors  using  two-level  caches,"  CSL  Kept. 
CSG-36,  Univ.  of  I1L,  1984.  (JSEP) 

2.  A.  Ram  and  JJL  Patel,  "Parallel  garbage  collection  without  synchronization  overhead,"  CSL  Rept. 
CSG-35,  Univ.  of  I1L,  1984.  (JSEP) 


Theses: 

1.  D.  Colglazier,  "A  performance  analysis  of  multiprocessors  using  two-level  caches;"  MS.  thesis;  Elec, 
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SCIENTIFIC  OBJECTIVES: 

This  unit  will  focus  on  the  fundamental  issues  in  the  design  and  analysis  of  fault-tolerant  com¬ 
puters.  As  computational  power  and  the  physical  complexity  of  computer  systems  increase,  so  does  the 
need  for  increased  reliability  and  fault  tolerance.  In  addition,  there  is  a  pressing  need  for  systems 
which  will  continue  to  function  under  a  variety  of  workload  environments  and  failure  conditions.  We 
have  identified  several  issues  which  require  particular  attention  due  to  both  technology  and  systems 
viewpoints.  These  approaches  are  necessary  for  realizing  effective  fault  tolerance  without  extreme  cost 
and  performance  penalties. 

1.1  Hierarchical  testing  and  algorithm-based  fault  tolerance 

Many  computer  systems  are  being  designed  with  multiple  units  (processors  and  memories)  which 
cooperate  to  perform  various  tasks.  Each  of  the  units  is  usually  quite  complex  and  is  designed  in  a 
hierarchical  fashion.  We  propose  to  exploit  these  characteristics  of  computer  systems  in  deriving  tech¬ 
niques  to  analyze  and  design  fault-tolerant  systems. 

We  propose  to  derive  techniques  which  will  enable  the  results  obtained  at  any  level  in  a  system 
to  be  applied  to  higher  levels  in  the  hierarchy.  Currently,  no  useful  techniques  exist  which  are 
hierarchically  decomposable.  Thus,  for  example,  many  techniques  exist  to  derive  test  patterns  for 
modules.  However,  given  a  collection  of  modules  (with  associated  test  patterns)  and  their  interconnec¬ 
tion,  it  is  not  easy  to  derive  tests  for  the  collection.  The  same  problem  exists  in  applying  fault-tolerance 
techniques  such  as  error  detection  or  recovery  to  higher  levels. 

Another  objective  of  this  research  is  to  study  fault-tolerant  algorithms  for  systems  consisting  of 
multiple  intelligent  units.  Editing  techniques  for  fault  tolerance  are  largely  brute  force,  involving 
masnve  amounts  of  redundancy  in  schemes  such  as  triplication  and  voting.  Utilizing  the  knowledge  of 
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the  algorithms  which  execute  on  the  system  as  well  as  the  knowledge  of  the  intelligent  units  in  the 
system  and  their  interconnection  can  result  in  fault  tolerance  techniques  which  require  only  negligible 
overhead  in  hardware  or  performance.  This  novel  approach  to  fault  tolerance  will  be  studied  in  detail 
in  this  unit. 

We  will  also  develop  a  formal,  theoretical  foundation  for  the  analysis  and  design  of  fault- 
tolerant  systems.  This  theory  will  allow  the  designer  to  decide  just  how  much  and  what  form  of 
redundancy  is  appropriate  for  a  given  functional  module  in  a  given  environment.  The  analysis  tools 
will  enable  the  study  of  trade-offs  between  performance  and  reliability,  and  the  relation  of  these  to 
system  complexity. 

1.2  Measurement-  and  symptom-based  analysis  of  f aidx-toleranx  systems 

Modeling  hardware  and  software  fault  tolerance  is  a  complex  analytical  problem.  Results  based 
on  real  measurements  and  experiments  are  essential  if  a  scientific  basis  is  to  be  developed  for  accurate 
analytical  modeling.  This  research  is  concerned  with  the  development  of  effective  reliability  models  to 
describe  both  the  performance  and  reliability  aspects  of  computer  systems.  This  work  is  primarily 
motivated  by  the  fact  that  previous  research  shows  that  system  reliability  is  a  dynamic  function  of  the 
level  and  nature  of  system  activity  as  described  by  the  system  workload.  These  results  have  major 
significance  because  they  indicate  that  (in  contrast  with  conventional  wisdom)  it  is  not  useful  to  push  a 
system  close  to  its  performance  limits;  since  the  gain  in  performance  is  more  than  offset  by  degradation 
in  reliability.  The  proposed  research  is  in  two  related  directions. 

First,  we  propose  to  develop  hardware  and  software  reliability  models  to  explicitly  take  the 
nature  and  level  of  system  activity  into  account.  These  models  will  be  based  on  real  error  and  work¬ 
load  data  being  collected  on  several  systems  on  campus.  IBM,  DEC  and  Cyber  systems  are  being  moni¬ 
tored  to  provide  appropriate  data  for  both  model  development  and  validation. 

Second,  it  is  proposed  to  study  possible  cause-and-effect  relationships  to  model  the  observed 
workload/failure  dependency.  Two  important  investigations;  in  this  context,  are:  (a)  the  effect  of 
workload  variations  on  error/fault  latency  characteristics  and  (b)  the  reliability  impact  of  electrical 
and  thermal  stresses  induced  by  varying  levels  of  system  activity. 

We  also  propose  to  develop  theoretical  foundations  and  practical  models  for  inference-based  adap¬ 
tive  recognition  of  related  errors  in  large  computer  systems.  The  objective  is  to  investigate  the  feasibil¬ 
ity  of  recognizing  the  symptoms  of  severe  errors  based  on  relating  simpler  but  similar  observed  error 
symptoms.  The  validation  of  captured  error  symptoms  will  primarily  be  through  recurrence.  It  is 
envisaged  that  the  methodology  developed  will  allow  capturing  of  symptoms  (patterns)  of  severe  faults 
in  large  systems.  The  patterns  found  are  expected  to  form  the  basis  of  an  inference-based  error  predic¬ 
tion  methodology.  It  is  proposed  that  the  technique  developed  be  based  on  the  observed  error  rate  and 
the  type  of  activity  in  progress  at  the  time  of  the  error  and  validated  pattern  recurrence.  The  tech¬ 
nique  is  expected  to  be  general  so  as  to  encompass  different  machine  architectures  and  configurations. 

1.3  Reliable  algorithms  and  architectures  for  database  systems 

This  research  is  concerned  with  the  analysis  and  development  of  fault  tolerant  special  purpose 
architectures  and  algorithms  for  database  systems.  The  recent  interest  in  knowledge-based  fifth  genera¬ 
tion  computing  systems  has  resulted  in  new  challenges  concerning  special  purpose  database  system 
architectures.  Though  there  has  been  a  considerable  amount  of  research  concerning  large  scale  database 
systems  there  has  been  little  fundamental  work  on  architectures  and  algorithms  for  special  purpose 
fault-tolerant  multiprocessor  systems  appropriate  for  VLSI  implementation.  We  propose  to  derive  an 
integrated  approach  to  fault  tolerance  in  both  tightly-coupled  and  loosely  distributed  database  systems 
that  includes  novel  fault  tolerant  distributed  algorithms  and  reconfigurable  architectures  with  high 
performance  reliable  operation  in  harsh  environments.  Previous  approaches  to  fault  tolerance  in  such 


systems  have  relied  on  system  crashes  or  vague  notions  of  assertion  checking  to  detect  errors  during  sys¬ 
tem  operation.  Our  approach  is  to  use  algorithm-specific  architectural  enhancements  to  detect  module- 
level  failures  and  allow  real-time  system  level  reconfiguration  and  recovery. 


SUMMARY  OF  RESEARCH: 

2.1  Algorithm-based  fault  tolerance 

The  concept  of  algorithm-based  fault  tolerance  deals  with  system-level  methods  for  obtaining 
reliable  results  from  computations,  especially  when  performed  on  array  processor  systems.  In  this 
scheme,  algorithms  have  their  outputs  encoded  in  a  system-level  error-detecting  or  -correcting  code  and 
provide  appropriate  error-checking  procedures.  We  have  performed  a  probabilistic  analysis  of  the 
scheme.  Expressions  for  the  reliability  of  the  results  of  the  computation  and  the  time  for  completion  of 
the  computation,  using  a  particular  algorithm,  have  been  derived  in  terms  of  various  parameters:  the 
number  of  processors  involved,  the  time  for  execution,  and  the  fault-detecting,  -locating,  and  -tolerating 
capabilities  of  the  algorithm.  Using  a  semi-Markov  process,  we  have  modeled  the  system  with  con¬ 
current  error  detection,  with  recovery  and  reconfiguration.  The  model  can  handle  any  failure  distribu¬ 
tion  of  the  processors  and  any  reconfiguration  and  retry  time  distribution.  The  only  restriction  is  that 
it  is  assumed  that  the  time  for  executing  a  single  instance  of  a  problem  is  fixed.  This  simplification  can 
be  justified  if  the  algorithms  have  the  property  that  the  execution  time  depends  on  the  problem  size  and 
not  the  values  of  the  data  elements.  Hence,  given  a  certain  problem  size,  the  algorithm  always  takes  a 
fixed  amount  of  time.  For  mast  array  processor  algorithms  which  are  characterized  by  a  relatively  sim¬ 
ple  control  flow,  this  assumption  is  valid.  This  measure  can  aid  in  comparing  various  alternate  algo¬ 
rithms  to  solve  a  given  problem  if  the  selection  is  made  on  the  basis  of  the  reliability  of  the  produced 
results.  Another  measure  of  effectiveness  of  a  fault-tolerance  scheme  that  has  been  proposed  is  related 
to  the  number  of  problems  that  can  be  solved  on  a  system  with  a  bound  on  the  number  of  processors 
that  can  be  used  and  a  bound  on  the  reliability  of  the  result.  This  measure  is  restricted,  however,  to 
problems  that  can  be  partitioned  into  smaller  subproblems. 

2.2  Measurement-  and  symptom-based  analysis  of  fault-tolerant  systems 

Our  preliminary  investigations  have  resulted  in  the  development  of  a  general  methodology  to 
determine  the  effect  of  workload  variations  on  error  latency.  This  technique  was  applied  to  determine 
the  effect  of  latent  memory  errors  on  a  VAX  11/780  system.  Preliminary  results  on  thermal  stresses 
due  to  system  activity  show  a  strong  contribution  due  to  supervisor  state  activity. 

We  have  implemented  a  methodology  to  detect  and  relate  automatically  the  pattern  for  frequent 
errors  on  two  Cyber  systems.  This  technique  captures  symptoms  of  severe  faults,  including  those 
which  propagate  across  coupled  machines. 

2.3  Reliable  algorithms  and  architectures  for  database  systems 

Our  initial  investigations  have  resulted  in  two  novel  methods  of  incorporating  redundancy  in 
linked  data  structures  which  allow  concurrent  checking  of  structural  integrity  during  normal  access 
path  traversal.  Our  first  proposal  is  referred  to  as  distributed  and  appended  checks  and  is  based  on 
code-theoretic  concepts.  The  second  proposal  incorporates  compressed  structural  checks  in  the  form  of 
signatures  and  is  called  signatured  access  path  checking.  In  contrast  to  the  few  suggestions  made  by 
previous  researchers,  our  results  provide  a  unified  approach  to  real-time  error  detection  in  large  database 
systems. 
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SCIENTIFIC  OBJECTIVE: 

Current  technologies  require  the  development  of  efficient  computational  techniques  and  the 
analysis  of  the  capabilities  of  various  models  of  computation.  We  are  concerned  with  the  resources  - 
such  as  time,  equipment,  memory,  interconnection  -  that  are  used  or  needed  in  the  algorithmic  solution 
of  given  problems.  This  dynamic  discipline,  known  as  concrete  computational  complexity,  not  only 
contributes  to  a  fundamental  understanding  of  computing,  but  is  also  extremely  relevant  to  both 
hardware  and  software  applications.  It  is  only  natural  that  the  great  technological  innovations 
represented  by  Very-Large-Scale-Integration  (VLSI)  and  the  advent  of  computer  networks  have  had  a 
substantial  impact  on  this  field,  opening  new  horizons  and  posing  challenging  problems.  Much  of  our 
current  research  (organized  below  in  six  subsections)  has  been  motivated  by  these  innovations. 
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SUMMARY  OF  RESEARCH: 

1.  Computational  aspects  of  VLSI 

Our  research  in  this  area  can  be  broadly  subdivided  into  two  main  categories:  VLSI  Algorithms 
and  Layout  Theory. 

12  VLSI  Algorithms 

A  dictionary  is  a  data  structure  that  supports  insertion,  deletion,  and  retrieval  operations.  To 
maintain  a  database,  a  dictionary  machine  accepts  an  arbitrary  sequence  of  instructions  at  a  constant 
rate.  We  designed  new  VLSI  dictionary  machines  on  a  shuffle-exchange  network  (SEN)  and  on  a  cube- 
connected  cycles  network  (CCC).  The  SEN  design  has  a  novel  pipelined  architecture.  Each  stage  of  the 
pipeline  operates  on  batches  of  instructions,  not  on  individual  instructions  The  SEN  runs  forward  to 
execute  instructions,  backward  to  deliver  the  responses  to  the  instructions.  To  implement  a  dictionary 
machine  on  a  CCC  we  gave  the  first  explicit  proof  that  every  CCC  has  a  Hamiltonian  path. 

With  the  intent  to  take  a  further  step  toward  the  development  of  a  coherent  theory  of  VLSI  com¬ 
plexity,  we  have  concentrated  on  the  information  exchange  that  takes  place  in  a  chip  during  the  execu¬ 
tion  of  an  algorithm.  In  particular,  we  have  investigated  subdivisions  of  the  chip  that  are  more 
demanding  on  the  information  flow  than  the  traditional  balanced-bisection  technique.  The  novel  tech¬ 
nique,  called  "square  tessellation",  partitions  the  chip  into  a  collection  of  square  cells,  subsumes  bisection, 
and  in  many  cases  leads  to  stronger  lower  bounds  on  the  AT2  measure;  moreover,  it  reveals,  in  chips 
with  a  high  input  density  (measured  in  bit/area),  a  new  phenomenon  —  called  "saturation"  — ,  where 
storage  limitations  pose  stronger  requirements  on  the  information  flow.  We  have  obtained  new,  more 
powerful  general  lower-bound  techniques,  which  we  have  successfully  applied  to  problems  such  as 
sorting  and  shifting.  In  addition  we  have  explored  the  complementary  upper-bound  problem  and  pro¬ 
duced  a  full  spectrum  of  optimal  or  near-optimal  parallel  architectures  for  the  above  problems. 

13  Layout  Theory 

Under  the  Thompson  model,  it  is  well-known  that  bounded-degree  n-node  graphs  with  0('/n) 
separator  can  be  laid  out  in  linear  area  with  maximum  edgelength  (X>/n/lognX  and  n-node  bounded- 
degree  graphs  with  F=fi(>/n)  bif urea  tor  can  be  laid  out  in  CXF2log2Cn/F))  area  with  maximum 
edgelength  0(Flog(n/F)/loglog(n/F)).  These  results  are  also  existentially  optimal  (to  within  a  constant 
factor)  in  the  sense  that,  for  each  class,  there  exists  a  graph  in  the  class  that  requires  the  specified 
amount  of  area,  and  has  the  specified  maximum  edgelength. 

We  generalize  and  unify  these  results  by  considering  embeddings  in  variable  aspect-ratio  bound¬ 
ing  rectangles.  We  show  that  if  we  consider  minimum-area  embeddings  of  classes  of  graphs  with 
separator  F=o(  Vn)  or  bifurcator  F*fl(VnX  then,  if  the  longer  side  of  the  bounding  rectangle  is  a  and 
the  shorter  side  is  b,  there  exists  a  graph  in  the  class  for  which  the  maximum  edgelength  in  the  embed- 
ing  is  Q(a/log(b/F)X  For  the  special  case  when  (squareroot  of  the  area),  this  gives  the  previ¬ 

ously  known  results.  Further,  our  proof  is  simpler  than  the  earlier  proof  for  the  lower  bound  on  the 
maximum  edgelength  for  bif  urea  tor-based  embeddings,  and  the  results  are  existentially  optimal,  since 
we  have  exhibited  layout  strategies  that  achieve  this  bound  on  the  maximum  edgelength. 

Frequently  in  the  context  of  layout  theory  one  must  face  a  variety  of  geometric  retrieval  prob¬ 
lems.  We  have  developed  an  efficient,  general  search  technique  (called  "groping”  search)  that  applies  to 
problems  for  which  only  inefficient  solutions  were  known,  such  u  circular  and  half-space  retrieval 

2.  Distributed  Computation 

We  have  continued  to  study  fundamental  problems  in  distributed  computation. 

First,  we  devised  an  election  algorithm  on  a  complete  network  of  n  processors  that  has  a  sense  of 
direction  with  Ofn)  messages.  In  contrast,  it  is  known  that  on  a  complete  network  without  a  sense  of 
direction  the  masssge  complexity  of  election  is  Ofn  log  n). 
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Second,  we  characterized  completely  the  shared  memory  requirements  for  achieving  agreement  in 
an  asynchronous  system  of  n  fail-stop  processes  that  die  undetectably.  There  is  no  agreement  protocol 
that  uses  only  read  and  write  operations.  This  result  implies  the  impossibility  of  Byzantine  agreement 
in  asynchronous  message-passing  systems,  which  had  been  established  by  MJ.  Fischer,  N.A.  Lynch,  and 
M.  Paterson  under  slightly  different  conditions.  Furthermore,  there  is  an  agreement  protocol  with  test- 
and-set  operations  if  and  only  if  either  memory  cells  have  at  least  three  values  or  n  »  2. 


INTERACTION  WITH  OTHER  RESEARCHERS: 

At  Illinois  we  expect  to  work,  closely  with  E  E  Hajek  in  CSL,  D.  E  West  in  the  Department  of 
Mathematics,  and  C  L.  Liu  in  Computer  Science.  Professor  Hajek’s  activity  in  computer  communica¬ 
tion  networks  provides  a  source  of  relevant  problems  for  distributed  computation.  Professors  West  and 
Liu  provide  expertise  on  the  combinatorial  aspects  of  our  work.  Our  research  in  all  of  the  above  areas  is 
related  and  complements  analogous  activity  in  major  research  institutions;  indeed,  we  communicate  reg¬ 
ularly  with  researchers  at  MJLT,  Yale  University,  Berkeley,  Cornell  University,  Carleton  University, 
University  of  Rochester,  and  IBM  Research  /  San  Jose. 
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SCIENTIFIC  OBJECTIVE: 

The  term  multi-sensor  array  system  refers  to  a  large  class  of  remote  sensing  systems  in  which  data 
are  collected  and  recorded  by  many  independent  sensors,  or  by  one  sensor  that  is  moved  to  different  spa¬ 
tial  positions.  The  recorded  data  are  processed  by  a  digital  array  algorithm  to  produce  a  high  resolution 
object  function.  Some  of  the  more  important  multi-sensor  array  systems  now  in  use  are  synthetic  aper¬ 
ture  radar  (SAR),  computer-aided  tomography  (CAT)  and  beamforming  sonar.  The  objective  of  the 
research  in  this  unit  is  to  develop  both  the  theory  and  computer  verification  of  new  signal  processing 
methods  with  the  goal  of  overcoming  current  limits  on  resolution  and  data  throughput  rates  for  these 
systems. 


SUMMARY  OF  RESEARCH: 

During  the  past  year  our  JSEP  research  has  concentrated  on  band-limited  signal  extrapolation  and 
phase  retrieval,  multidimensional  windowing,  and  on  the  development  of  the  convolution  back- 
projection  image  reconstruction  algorithm  for  spotlight  mode  SAR. 


SIGNIFICANT  RESEARCH  ACCOMPLISHMENTS: 

Both  band-limited  signal  extrapolation  and  phase  retrieval  are  techniques  for  reconstructing  sen¬ 
sor  array  images  from  incomplete  observations.  Our  studies  of  these  techniques  have  produced  a 
number  of  major  results,  including: 


(1)  Alternative  formulation  of  the  extrapolation  problem  in  the  presence  of  noise,  which  may  lead 
to  robust  algorithms. 

(2)  Linking  results  on  extrapolation  in  the  signal  processing  literature  to  those  in  mathematics 
(integral  equations,  Hilbert  space). 

(3)  Two  generalized  discrete  band-limited  extrapolation  algorithms  which  are  much  more  flexible 
than  existing  algorithms. 

(4)  Characterization  of  algorithms  for  band-limited  extrapolation  using  a  linear  time-varying  model 
that  clearly  shows  why  certain  algorithms  work  well  for  some  examples  and  fail  for  others. 

(5)  Proving  that  the  phase  retrieval  problem  is  both  stable  and  ill-conditioned. 

In  many  multi-sensor  problems  only  a  finite  region  of  the  spatial  or  frequency  domain  can  be 
sampled.  In  previous  work  we  developed  a  computer  optimization  technique  for  generating  "optimal" 
windows  that  can  be  used  to  improve  image  quality  when  finite  areas  of  irregular  shape  are  processed 
[lj.  Windows  designed  by  this  technique  are  optimal  in  the  sense  that  as  much  energy  as  possible  is 
concentrated  within  a  narrow  mainlobe  in  the  spectral  domain,  while  maintaining  the  finite  and  irreg¬ 
ularly  shaped  support  in  the  spatial  domain.  While  this  procedure  produces  useful  windows  that  have 
been  shown  to  improve  image  quality,  the  design  procedure  is  numerical,  and  is  not  easily  characterized 
by  a  closed  form  analytic  expression. 

Our  work  on  this  problem  during  the  past  year  consisted  of  a  study  that  attempted  to  determine 
the  best  way  to  construct  good  2-D  windows  from  a  good  1-D  window  when  the  support  is  irregularly 
shaped.  The  approach  was  to  consider  a  small  set  of  2-D  geometric  shapes  (rectangles,  squares,  dia¬ 
monds,  circles,  L’s,  cones,  dual-cones,  etc.)  to  window  these  regions  by  placing  the  ID  window  along  cer¬ 
tain  selected  contours,  and  to  calculate  the  spectrum  of  the  result  by  means  of  a  2-D  EFT.  The  resulting 
patterns  were  characterized  by  the  mainlobe  widths  and  peak  sidelobe  responses  along  each  frequency 
axis,  and  along  the  45*  contours  in  the  frequency  plane.  After  analyzing  a  large  number  of  experimen¬ 
tal  windows  and  tabulating  the  results,  a  general  rule  suggesting  the  best  way  to  apply  the  1-D  win¬ 
dow  to  the  2-D  area  began  to  emerge.  The  rule  is  that,  to  the  extent  possible,  the  1-D  window  should  be 
placed  along  contours  which  meet  the  edges  of  the  support  at  right  angles.  For  rectangles,  the  1-D  win¬ 
dow  should  be  applied  along  perpendicular  lines  that  run  parallel  to  the  axes.  For  circles  the  1-D  win¬ 
dow  should  be  applied  along  each  diameter,  creating  a  circularly  symmetric  window.  On  an  annular 
sector,  the  1-D  window  should  lie  along  circular  contours  and  radial  contours,  etc.  Although  it  is 
difficult  to  form  an  analytical  basis  for  the  90  *  guideline,  this  rule  does  appear  to  be  useful  for  con¬ 
structing  2-D  windows  of  reasonable  quality  when  the  region  of  support  has  an  irregular  shape. 

In  some  of  our  past  work  we  were  able  to  show  that  there  is  a  similarity  between  computer-aided 
tomography  (CAT)  and  synthetic  aperture  radar  (SAR),  and  that  the  reconstruction  of  a  SAR  image  can 
be  interpreted  as  a  tomographic  reconstruction  based  on  a  coherent  limited-angle  recording  of  spatial 
domain  "projections"  or  Fourier  domain  "slices"  [2J.  This  discovery  led  to  the  suggestion  that  convolu¬ 
tion  back-projection  (CBP)  reconstruction,  which  is  a  state-of-the-art  algorithm  in  medical  imaging, 
could  also  be  used  for  SAR.  The  advantages  of  using  CBP  for  SAR  are  that  the  algorithm  operates 
directly  on  polar  data,  thereby  eliminating  the  need  for  error-prone  polar-to-cartesian  interpolation  [3l 
Also,  CBP  is  known  in  CAT  for  being  relatively  insensitive  to  noise  in  the  recording.  The  disadvantage 
is  that  CBP  reconstruction  is  a  point-by-point  reconstruction  algorithm  that  has  an  order  of  complexity 
of  O  (N^X  where  N  x  N  is  the  size  of  the  final  image.  In  contrast,  the  order  of  complexity  for  FFT- 
based  Fourier  domain  algorithms  is  0(kN2),  where  k  is  a  function  of  the  order  of  the  polar-to-cartesian 
interpolators.  Since  k  <  <N  in  typical  problems,  the  CBP  algorithm  requires  more  computation  than 
the  FFT  algorithms. 

However,  several  results  have  been  obtained  which  suggest  that  CBP  may  be  a  viable  SAR  algo¬ 
rithm  having  a  parallel  structure  which  is  more  suitable  for  VLSI  implementation  than  the  batch  FFT 
algorithms  First,  it  has  been  shown  by  computer  experiments  for  point  target  responses  that  the  CBP 
algorithm  produces  high  quality  reconstructions,  and  in  some  cases  eliminates  the  production  of  spuri¬ 
ous  responses  that  originate  in  the  polar-to-cartesian  interpolators.  Also,  it  has  been  discovered  that  the 


CBP  algorithm  can  be  more  easily  implemented  in  real-time  and  in  step  with  the  data  gathering  process. 
This  is  because  it  can  make  use  of  the  interpulse  transmission  intervals  for  the  filtering  of  the  projec¬ 
tions  (convolution),  and  can  carry  out  the  back-projection  in  real-time  as  the  recording  evolves.  Also, 
because  the  CBP  algorithm  is  a  point-by-point  algorithm  (rather  than  a  batch  algorithm)  it  can  be 
implemented  by  many  small  processors  working  in  parallel  to  produce  a  sub-area  of  the  final  image. 
This  suggests  a  type  of  systolic  array  architecture  which  is  attractive  for  VLSI  implementation.  Also,  a 
multiprocessor  architecture  has  been  proposed  that  offers  a  degree  of  fault  tolerance  that  can  be 
achieved  with  the  inclusion  of  error  detection  and  the  proper  management  of  the  total  resources. 
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SCIENTIFIC  OBJECTIVE: 

The  objective  of  the  research  is  to  tackle  fundamental  problems  of  electronic  materials,  devices 
and  systems  in  a  timely  manner,  and  to  provide  early  funding  on  start-up  of  projects  which  present 
immediate  opportunities  of  high  scientific  merit. 

SUMMARY  OF  RESEARCH: 

The  funds  in  this  unit  were  used  to  initiate  a  new  project  on  the  laser  growth  of  films  under  the 
direction  of  Professor  J.G.  Eden.  Graduate  student  GC.  Abele  has  been  studying  the  effect  of  illuminat¬ 
ing  surfaces  with  visible  or  UV  laser  radiation  during  the  CVD  growth  of  films.  A  more  detailed  sum¬ 
mary  of  this  research  appears  in  Unit  11  under  the  subtitle  "Laser  Growth  of  Films." 

In  addition,  these  funds  are  supporting  one-half  the  salary  of  a  research  associate,  Scott  Barnett. 
Scott  Barnett  recently  received  his  Phi),  degree  under  the  direction  of  Professor  J.E.  Greene.  Dr.  Barnett 
agreed  to  accept  a  post-doctorate  position  until  next  fall  in  order  to  maintain  continuity  in  their 
research  program  and  to  complete  several  unfinished  projects.  More  details  of  their  work  appears  in 
Unit  5. 


\  *,  s 


-  '-T5-”-  A-'-*-'.  "-'■'--•'-I-'- 


END 


FILMED 


7-85 


DTIC 


*  •^*r.  *" '* .  ^ ■*,  **,  •*,  »*,  «*, 

*  •/  •, **  %'"»/•  *  «  **».V 

•  *  -  •  *  •  »  *  #  *  U  *  J»  *  M  *  *  *  .  *  m  *  _  •  .  «  4  •  _  . 


